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■ Abstract 

^~>' In this paper novel classes of 2-D vector-valued spatial domain wavelets are defined, and their properties given. 

The wavelets are 2-D generalizations of 1-D analytic wavelets, developed from the Generalized Cauchy-Riemann 
equations and represented as quaternionic functions. Higher dimensionality complicates the issue of analyticity, more 
[ than one 'analytic' extension of a real function is possible, and an 'analytic' analysis wavelet will not necessarily 

ps| ■ construct 'analytic' decomposition coefficients. The decomposition of locally unidirectional and/or separable variation 

is investigated in detail, and two distinct families of hyperanalytic wavelet coefficients are inttoduced, the monogenic 
and the hypercomplex wavelet coefficients. The recasting of the analysis in a different frame of reference and its effect 
on the constructed coefficients is investigated, important issues for sampled transform coefficients. The magnitudes of 
the coefficients are shown to exhibit stability with respect to shifts in phase. Hyperanalytic 2-D wavelet coefficients 
enable the retrieval of a phase-and-magnitude description of an image in phase space, similarly to the description 

■ of a 1-D signal with the use of 1-D analytic wavelets, especially appropriate for oscillatory signals. Existing 2- 
D directional wavelet decompositions are related to the newly developed framework, and new classes of mother 
wavelets are introduced. 
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O ri~l HIS paper constructs new classes of 2-D vector-valued mother wavelet functions. The vector-valued wavelet 
^ M. functions are defined to produce 2-D extensions of 1-D analytic complex wavelet coefficients, denoted hyper- 
analytic wavelet coefficients. The extension of 1-D analytic wavelet coefficients to 2-D is non-trivial. Firstly more 
^ , than one definition of 'analytic' exists in 2-D. The suitable choice of the concept in 2-D depends on the local 
' directionality and dimensionality of the signal under analysis, and local unidirectionality, as well as separability, 
will explicitly be considered. Secondly, unlike the 1-D case, a distinction can be made between wavelet coefficients 
obtained from using wavelets that are 'analytic', wavelet coefficients that are 'analytic' in their spatial index, and 
wavelet coefficients of 'analytic' 2-D functions. In 2-D it will furthermore become necessary to introduce the notion 
of ^-'analytic', or 'analytic' in a rotated frame of reference, to interpret the properties of the wavelet coefficients. 
We shall demonstrate that it is possible to choose a definition of 'analytic' suitably so that the hyperanalytic wavelet 
coefficients have tractable and useful properties in the rotated frame. The mathematical framework developed for 
the construction of the hyperanalytic wavelet coefficients also provides a context to some of the most promising 
new complex/vector- valued wavelet filter decompositions in the literature [1], [2], as well as allows us to construct 
new classes of vector-valued wavelets. 

Over twenty year ago, Grossman and Morlet [3] developed the Continuous Wavelet Transform (CWT) [4], using 
continuous complex-valued mother wavelets. Initial analysis based on wavelet decompositions was implemented 
using such mother wavelets. Both magnitude and phase descriptions of non-stationary signals were determined, and 
an early example of analysis include wavelet ridge methods proposed by Delprat et al. [5]. However subsequently 
for many years interest focused on the Discrete Wavelet Transform (DWT) and signal estimation. The DWT was 
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developed to implement the wavelet transform of time-compact mother wavelets and as compact discrete wavelet 
filters cannot be exactly analytic [6], real wavelets were used. A revival of interest in later years has occurred in 
both signal processing and statistics for the usage of complex wavelets [1], [7], and in particular complex analytic 
wavelets [8]-[ll]. This revival of interest may be linked to the development of complex- valued discrete wavelet 
filters [12] and the clever dual filter bank [8]. The complex wavelet transform has been shown to provide a powerful 
tool in signal and image analysis [1], where most of the properties of the transform follow from the analyticity of 
the wavelet function. This paper strives to build upon the strengths of recent developments, by deriving large classes 
of wavelets generalizing the concept of a 1-D local complex- valued analytic decomposition to a 2-D vector- valued 
hyperanalytic decomposition. 

The wavelet decomposition is in 2-D only a single example of a local decomposition. Local decompositions of 
an image in terms of spatial structure associated with a given scale and spatial position have with great success 
formed the basis for many different procedures of local image analysis and estimation [4], [13], [14]. As the 
variational structure in 2-D is much richer than in 1-D, the term 'local' may denote a variety of different forms 
of spatially limited structures, and so examples of decompositions include i.e. wavelets [4], curvelets [15] and 
bandelets [16] etc. A local decomposition method should be chosen with easy interpretability of the decomposition 
coefficients and compression in mind, as both of these characteristics facilitate the analysis of observed images. 
For local image analysis, the local directionality and dimensionality of variation are key features of the local visual 
appearance of the image. These two features need to be characterised, and can be used for further processing, 
such as image segmentation, feature detection and disparity estimation, see for example work by von Billow and 
Sommer [17], as well as that of von Btilow [18], Hwang et al. [19], Antoine et al. [20], [21] etc. To represent both 
local dimensionality and directionality at each spatial point more than a scalar-valued representation is required. 
Based on a fixed local dimensionality the characterisation of the image will of course vary. For example, images 
that are locally 1-D, are well characterised by the direction of the image variation, its magnitude and the scale of 
variation in the given direction. Texture features present in an image may correspond to locally separable structures, 
characterised by the direction of the separable structure, the local magnitude of variation and the scale of variation 
in two directions. Single texture components often correspond to variations associated with the same scale in both 
perpendicular directions. 

The observational axis of an image cannot be assumed to coincide with the natural axis of the local structure. It is 

therefore necessary that the local representations should not substantially change form with rotations and translations 
of the spatial index, but exhibit equivariance to such operations. Of course the CWT exhibits equivariance when 
calculated in continuous space at all values of the locality index [13], but if the transformation is to be discretized, 
it needs to satisfy some further constraints so that the representation is stable. In particular we require that the 
magnitude of the vector-valued decomposition should be stable to shifts in phase of the signal, or a sampled 
collection of coefficients of the process may exhibit high translation variance in comparison to a sampled collection 
of coefficients obtained from the same image subjected to a small spatial shift in both axes simultaneously [1], [22]. 
Translation variance is clearly an undesirable feature of a decomposition. The simple characterisation of any shift 
in index is also important for analysis of multiple images. Vector-valued representations, as will be demonstrated, 
can with a suitable choice of decomposition filter, both provide decompositions stable to small spatial shifts (or 
phase shifts) and allow for the easy characterisation of rotations and translations of the spatial index. 

As examples of successful existing vector-valued local representations of image structure we note the 2-D analytic 
signal [23], or the complex CWT [1], [22]. The 1-D analytic signal at a given value of xi provides a meaningful 
magnitude and phase description of a signal, if only a single oscillatory component is present. For each fixed 
scale the complex wavelet coefficients correspond to spatially local representations of the image in terms of two 
wavelet coefficients, but the magnitude and phase of the complex wavelet coefficients are not when an arbitrary 
complex-valued wavelet is used either meaningful nor are they interpre table. For this reason complex wavelet 
transform coefficients obeying suitable conditions are generally used [1], [24]. We aim to generalise the usage of 
spatially local two-vector representations to local three- and four-vector representations, using hyperanalytic wavelet 
decompositions. 

Vector-valued representations of local structure in terms of quatemionic representations, or equivalently three- 
and four-vectors, have been developed by von Biilow and Sommer [17], Felsberg and Sommer [25] as well as Hahn 
and Snopek [26]. In particular Hahn [27], von Biilow and Sommer [17], as well as Felsberg and Sommer [25], 
have developed 2-D extensions of the 1-D analytic signal. These representations give descriptions of local variation 
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in terms of a local magnitude and phase/s. Phase should represent local variations and structure, whilst magnitude 
is usually interpreted as local signal presence [28]. A local well-constructed magnitude and phase description of 
a signal is often very informative. However if more than a single narrowband component is present at a given 
spatial location in an image, the 2-D 'analytic' signal, no matter which of the above generalizations are chosen, 
will completely lack interpretation, just hke in the case of the 1-D analytic signal of a multi-component signal [29], 
[30]. 

The work by Hahn and Snopek [26] consists of forming a suitable 2-D quatemionic (four- vector valued) extension 
of the Wigner-Ville distribution [30]. In 1-D the Wigner-Ville distribution is known to be uninformative when 
describing signals contaminated by noise [31], and is unsuitable as a description of signals consisting of more 
than a single component at any spatial position due to interference between components [30]. To deal both with 
the addition of noise and the presence of multiple components in 1-D, local analytic representations have been 
constructed, such as analytic wavelet coefficients [5], [32], [33], or other local representations of the signal calculated 
in conjunction with a localisation procedure [34], [35]. In this paper we aim to construct 2-D 'analytic wavelet 
coefficients' or, locally 'analytic' representations of images, by using vector-valued wavelet functions. The wavelets 
will be constructed to merge the 2-D wavelet transform with existing 2-D 'analytic' representations. The synthesis 
of local image representations and magnitude and phase descriptions of local structure will be shown to enjoy many 
special properties, but requires several non-trivial generalizations of existing 1-D concepts. 

Before proceeding to the construction of 2-D wavelets to serve the same purpose as 1-D analytic wavelets, let 
us first briefly revisit the properties of the 1-D analytic wavelet decomposition. The decomposition is constructed 
from a single mother wavelet function ijj{xi) G (M) , that satisfies the admissibihty condition [13]. The mother 
wavelet is translated in time by h, and scaled by factor a to form a family of wavelets denoted by {i^a.bi^i)} ■ For 
V'(x) G (M*^) , the dilation and scaling operators for any index dimension d are Vai^ (x) = jap'^^^ (x/a) , and 
7b (x) = V' (x — b) , respectively. A member of the 1-D wavelet family is ipa,b{xi) = T^a^i^ (xi) ■ To represent 
a 1-D signal g{x) in terms of local contributions the CWT, and its inversion, are defined by 

POO l-OD l-OO 

w^{a,b;g)= c/(a;i)V'*,b(xi) dxi, g{xi) = — / w^{a,b; g)ipa,b{xi) —db. (1) 

J —oo ^tp J —oo J —oo ^ 

e M is a constant specific to the function ^^{xi), its existence guaranteed by the admissibihty condition. If '0(-) 
can be interpreted as local in time to xi = and frequency to /i = /o, the signal is reconstructed in terms of the 
local contributions w.^{a,b;g) : furthermore w^{a,b;g) is associated with a function local to position xi = b and 
frequency fi = fo/a. The wavelet coefficients are used to characterise the local properties of the signal and the 
decomposition is often described as a 'local Fourier' transform, even if this is a somewhat inaccurate description. 

If ip{xi) is real-valued, then the magnitude of w^{a,b]g) can be used to determine local signal content, or be 
used as the basis of more complicated analysis methods. If w^{a, b; g) is complex-valued then it can be represented 
by its magnitude and phase. If '0(a;i) is analytic, denoting an analytic mother wavelet ^ by il)^'^\xi), then the 
wavelet is a complex-valued function. The CWT of a real-valued signal using a complex-valued analytic wavelet 
yields complex-valued wavelet coefficients, where the coefficients are an analytic signal in b. In 1-D, if we define 
the complex variable ti = xi + jyi and restrict yi > 0, then an analytic signal g^'^\xi) = g{xi) + jg^^\xi), 
evaluated in the complex argument, that is g^'^\ti), is an analytic function in the upper half-plane of yi > [23]. 
g^~^\ti) thus satisfies the Cauchy-Riemann equations in the upper half -plane and this determines the relationship 
between g{ti) and g^^\ti). The fact that a real and imaginary part of g^'^\xi + jyi) satisfy the Cauchy-Riemann 
equations imply that the phase of g'^'^\xi) will be well defined, once some additional restrictions are placed on 
g{xi). Given the analytic wavelet coefficients are an analytic signal, and local in frequency, the instantaneous 
frequency of an oscillatory component present at scale a, may for example be determined using its analytic wavelet 
coefficients, and multi-component oscillatory signals may be characterised via wavelet ridge analysis [5]. Note 
that if w^{a,b;g) G C but V'(^ci) 7^ then w^{a,b;g) can still be represented in terms of a magnitude 

and phase, but the phase will not in general be interpretable in terms of local structure. Furthermore wavelet 
coefficients constructed from analytic mother wavelets have a magnitude invariant to phase-shifts and this alleviates 
observed problems with coefficients exhibiting strong shift-dependence [1], [36]. The analytic wavelet has one major 
drawback, longer essential time-support, but provides the framework for analysis shorter time-support filters may 

'An 'analytic' mother wavelet in this context is a complex-valued mother wavelet that is also an analytic signal [30]. 
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approximately provide, see for example work by Selesnick et al. [1] as well as Femandes et al. [24]. Selesnick et al. 

give a comprehensive discussion of the many important properties of analytic wavelet coefficients when calculated 
in both 1-D and higher dimensions, and these properties justify our interest in developing a single mathematical 
framework for the interpretation of 2-D locally 'analytic' representations. 

In 2-D the CWT of (/(x) € L^(M^) is implemented by choosing an admissible mother [13] wavelet '^(x) G 
L^(M^), and from this function, a family of functions {?/;^(x)} for ^ = [a,9,h]^ is defined by: 

^g(x) = VaReT^ m (x) = i7gV(x). (2) 

Reg^x.) = g{r^ox), is the rotation operator given in matrix notation as: 

re = [{cos(0), - sm{e)} , {sm{9), cos{9)}] . 

The purpose of including 9 in the decomposition is to identify local behaviour like unidirectional, or separable 
structures, not aligned with the observational axes by rotating the analysis function over the full set of possible 
orientation, when necessary [13]. The wavelet coefficients of function g{x) are then defined by [13]: 

/oo poo 
/ 5(x)V'|(x) (fx, 
-oo J —OD 

where * denotes the act of conjugation. Just like in 1-D the Wjp{^,g) coefficients can be used to reconstruct the 
function g{x) from a weighted average of the {V^^(x)} by: 

5(x) = — / / / / w^{^,g)^^{x) 3 . (3) 

J-oo J-oo Jo -'0 ^ 

To form Wjp{^,g), g{x) has been localised in both space and spatial frequency depending on the form of V'(x). 
^(x), the mother wavelet function, thus needs to be chosen with care, so that w^{^,g) has suitable properties. We 
wish to derive w^{^,g) that are locally 'analytic'. 

Firstly we need to specify what we mean by 'analytic' in 2-D. A reasonable starting point to generahse 
to 2-D is by considering 2-D versions of analytic functions and then forming the appropriate Umit, 
as limj^j_+o+ g^'^\h) = g^'^\xi). The definition of a 2-D analytic function has previously been examined in detail 
in pure mathematics, within the field of Harmonic Analysis, by Stein and Weiss, see for example [37], [38]. Similar 
concepts have appeared in applied mathematics and geophysics [39], [40], and in signal processing [17], [25], [27]. 
An important fact to note is that the generalisation of analyticity to higher dimensions is not unique [38], and that 
more than a single generalisation of the Cauchy-Riemann system, and thus to the definition of a hyperanalytic 
function, is possible. Any function that satisfies a generalization of the Cauchy-Riemann system will in fact be 
referred to as a hyperanalytic function in this paper. A hyperanalytic signal is, in contrast to the hyperanalytic 
function, defined as the suitable limit of such a function. 

Two possible hyperanalytic signals are relevant in this context: the hypercomplex signal of von Btilow and Sommer 
[17], and the monogenic signal of Felsberg and Sommer [25]. For clarity we will provide some discussion of the 
two hyperanalytic signals that are both possible and useful 2-D extensions of the 1-D analytic signal, corresponding 
in turn to four- and three-vector valued functions respectively. The hyperanalytic signal is usually represented as a 
quatemionic function, see for example [17], [25], i.e. as an object taking values in the skew field of quaternions, 
denoted H. The advantage of using the quaternion representation is that convenient polar forms (magnitude and 
phase/s) are naturally introduced [17], [25], representing magnitude and variation of either separable or unidirectional 
variation. Once the choice of 'analytic' has been specified, or convenient local polar representation, the 'analytic' 
mother wavelet is formed. We approach the choice of mother wavelet by striving to separate disparate structures 
present at the same spatial position and to represent the isolated local structure, in terms of local structure (phase) 
and magnitude. The interpretation of the phase functions of these coefficients must be established, and this will 
depend on the form of ip{x), and assumptions regarding the local structure of ^(x). We shall construct quatemionic 
wavelets starting from a real 2-D wavelet that will be augmented into a quatemionic object, in analogue to 1-D 
analytic wavelets. The constmction will rely on the choice of the starting point of a real mother wavelet, and the 
chosen form of hyperanalytic extension of this real function. 

One starting point to the constmction of the quatemionic mother wavelet is the anisotropic, or elongated, real- 
valued mother wavelet. Such a wavelet can constmct coefficients representing plane wave or separable stmctures 
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well. The wavelet would be extended into a quatemionic function so that a hyperanalytic wavelet coefficient 

representation can be formed of unidirectional or separable local structure in terms of the hyperanalytic phase/s 
and magnitude. Using quatemionic coefticients enables the easier parameterisation of local variation, and will rely 
on the transform coefficients being hyperanalytic, rather than the wavelet functions themselves. The strategy is 
related to work by Selesnick et al. [1] (who use directional complex wavelets), as well as Chan et al. [2] (who use 
separable quatemionic wavelets). Both the construction by Selesnick et al. and Chan et al. are implicitly based on 
the hypercomplex signal, and we explicitly construct a formal framework for interpreting such representations. We 
also form the monogenic extension of real directional wavelets. 

Another possibiUty is to separate the signal into components present at given scales using an isotropic real 
wavelet, and then assuming locally at that spatial position and scale either plane waves, or potentially separable 
structures only are present, represent these by forming the hyperanalytic signal of the wavelet coefficients. This 
corresponds to the 2-D version of the wavelet projection representations proposed by OUiede and Walden [34]. 

One of the main questions to tackle in the construction is the appropriate interpretation that can be given to the 
quatemionic wavelet coefficients depending on the choice of mother wavelet function. The answer is complicated 
in 2-D by the fact that rotations do not commute with the operators that construct hyperanalytic images. Given the 
complication introduced by the rotations, we start by considering analysis when no such operation is required: we 
assume that the image is aligned with the observational axes, and set the angle of rotation in the decomposition 
equal to zero. We then construct wavelet functions such that the CWT coefficients are hyperanalytic signals in 
their spatial indexing, and refer to such mother wavelets as hyperanalyticizing. The mother wavelet function that 
constracts monogenic coefficients, is a monogenic signal, but we refer to any quatemionic analysis filter that 
produces hypercomplex wavelet coefficients as a hypercomplexing wavelet, thus explicitly noting that the wavelet 
produces hypercomplex coefficients but is not necessarily itself a hypercomplex signals. This then, when no rotation 
is implemented, produces hypercomplex wavelet coefficients from a real-valued image. 

To deal with the act of rotation we introduce the definition of a ^-hyperanalytic signal, that is a function if 
observed in some rotated frame of reference, is a hyperanalytic signal. Once the monogenic signal has been 
extended to ^-monogenic, it follows that the quatemionic wavelet producing the ^-monogenic wavelet coefficients 
is itself ^-monogenic. In contrast, the hypercomplexing wavelet is not a 0-hyperanalytic signal as the property 
does not even hold for 9 = 0. We discuss in detail the interpretation of the wavelet coefficients constmcted 
from the family of wavelet functions constructed by scaling, rotating and translating the hyperanalyticizing mother 
wavelet. In analogue to the 1-D analytic and anti-analytic decomposition of a real- valued signal, see [41], we 
introduce the decomposition of a real- valued image into ^-hyperanalytic and anti-hyperanalytic components. The 
hyperanalyticizing wavelet annihilates/converts the anti-components to constmct an object equal to the real CWT of 
a 0-hyperanalytic function. This equivalence is important not only for the interpretation of the constracted coefficient 
but also allows for the derivation of additional useful properties of the hyperanalytic CWT. 

Specific features of local stmcture may be extracted directly and easily from the defined wavelet coefficients, for 
instance potential local unidirectionahty and separabiUty in the rotated frame of reference may be characterised. 
For aggregations of oscillatory signals, the wavelets we develop in this paper, can be used to separate the multiple 
components, and simultaneously give simple descriptions of the components. Previous work on local analysis of 
oscillatory signals is clearly distinct, a number of authors, see for example [19], [21], [42], [43], have discussed 
wavelet ridge analysis using complex wavelets. The advantage of using quatemionic rather than complex wavelets, 
as suggested in this paper, is that separable oscillations admit a simple and compact representation. Furthermore 
for local plane waves using a suitable choice of quatemionic wavelets, it is not necessary to calculate the transform 
at all directions to determine the local plane wave description. Thus the usage of quatemionic wavelet functions 
either introduces the simplified description of additional texture stmcture, or simplifies the analysis procedure. 

Coefficient equi- and invariance to local changes of phase and rotation alignment between the signal and the 
wavelet are also discussed. For the CWT this is not of immediate importance, as the transform coefficients 
exhibit equivariance with respect to changes in axes [13], but as discrete space decompositions may be related 
to continuous space decompositions, the stability of the decomposition to small spatial misalignment (phase-shifts) 
is very important [22]. For the hypercomplex and monogenic wavelet coefficients, we show that with a suitable 
choice of wavelet, magnitude invariance to phase shifts is achieved. For the monogenic wavelet coefficients, mother 
wavelet functions may be selected, such that the coefficients also exhibit invariance of magnitude with respect to 
local rotations. 



STATISTICS SECTION TECHNICAL REPORT: TR-06-02 



6 



We discuss generic structures of mother wavelet constructions, i.e. combining separable/ isotropic/ directional 
scale localising procedures, with both the hypercomplex and monogenic signal representation, but also give specific 
families of mother wavelets falling into these classes. Existing quaternionic filters can be related to the hypercomplex 
and monogenic constructions [2], [44]-[47]. In analogue to 1-D wavelets approximating a local Fourier transform, 
most of the existing quaternionic wavelets have been constructed to serve the role of a local Quaternion Fourier 
Transform (QFT). In addition to [2], work in a similar vein is due to [48]-[51] and corresponds to constructing 
local versions of the QFT. He & Yu [52], [53] constructed quaternionic multi-variate decompositions, and also 
work by Hsieh [54], [55] on motion smoothing should be noted, as well as work by Felsberg and Sommer on the 
monogenic scale space [47]. We provide a single complete framework that relates all of these constructions, and 
link the existing decompositions to Harmonic Analysis. We also construct completely new families of quaternionic 
wavelets that belong to the same family, and determine their properties. 

Central to the developments in this paper is the concept of hyperanalyticity and local versions thereof. The 
hyperanalytic signal is developed as a generalization of the 1-D analytic signal. The analytic signal is best introduced 
via the 1-D Fourier transform, 1-D local oscillations and analytic functions, and so we start by a discussion of 
these aspects of the analytic signal. 



II. Oscillations, dimensionality & analyticity 

A. Analyticity 

An oscillatory 1-D signal c(xi) is modelled as an amplitude and frequency modulated signal [29]: 

c(xi) = ac(a;i) cos(27r(^c(a::i)) = ^{adxi) exp{27rjipcixi))} . 

To determine the frequencies present in signal c(xi) the Fourier Transform (FT) of the signal could be calculated. 
The FT of a generic d-D signal ^(x), denoted G(f), is defined by: 

G(f) = / 5(x)exp(2j7rf^x) d^x, G(f) = |G(f)| exp(-2j7r(^G(f)), (4) 

and g{x.) can be reconstructed via 

g(x) = / |G(f)| exp(2j7r(f^x - (^G(f))) d'^i. (5) 

In 1-D at X = xi if g{xi) = c(xi) is oscillatory, and has a clearly defined content at frequency /i = fo{xi) 
approximately we then find c(xi) = |C(/o(xi))| cos(27r(/o(xi)j;i — ipc{fo{xi)))). To calculate the amplitude and 
phase of a real-valued signal, or |C(/o(xi))| and 27r(/o(xi)3;i — y'cC/o))) from c(xi), the analytic extension of a 
signal could be used [23]. This is defined for an arbitrary real signal g{xi) by: 



POD 

g^+\xi) = 2 / G(/)exp(27rjM) df = gixi) + jg^^\xi). 
Jo 



(6) 



g^^^ (xi) is the Hilbert Transform (HT) of g{xi). The HT in the time domain and the analytic signal in the frequency 
domain respectively, are 

g^^\x^)=n{g]{x^) = dy, gW(/i) = G (/i) (1 + sgn(/i)) . (7) 

J-oo xi-y 

The HT is as usual defined as a Cauchy principal value integral [23]. Using a stationary phase approximation [29] 
the analytic signal of c{xi) is then approximately c*^+)(xi) = c{xi) + jc^^\xi) = |G(/o(xi))| exp(2j7r(/o(xi)xi — 
^cifoixi)))) + o{\). Clearly the amplitude and phase of c{xi) are retrieved by adxi) = |c^^)(xi)| , and fdxi) = 
TT- tan^ ^1 (c(i)(xi)/c(xi)) . For a generic signal g{xi) the local ampUtude and phase may be defined by: 



ag{xi) 



g^^Hxi] 



2tt \ g{xi) 



For ag{xi) and ipg{xi) at point xi to be meaningful representations of g{xi), the latter signal is assumed to be at 
xi mainly limited in frequency to oscillations with a single period The key ingredient in defining a local amplitude 
and phase via equation ^ is the calculation of the analytic signal and the definition of g^^\xi). 
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The 1-D analytic signal is the limit of an analytic function, where a function is denoted analytic when it satisfies 
the Cauchy-Riemann equations in the upper half of the complex plane. To construct an analytic function in a complex 
argument from a real function in a real argument, the Poisson kernels and Poisson convolutions are introduced [38]. 
The Poisson kernels, and the convolutions of g{-) with these kernels, are using an auxiliary variable y defined by: 

1 y . 1 



PH{t 



1) = - 



Xl 



I) = - 



n xf + y'^ 



Ug{ti) = g{-) *PH{-,y), v^^\ti) = g{-) * qH{-,y)- 



(9) 



Define complex-valued variable ti = xi + jy, and note as y 0^, Ug{ti) —>■ g{xi) and vl^\ti) —>■ H{g} (xi). 



see [38] for details. kl^\ti) = Ug{ti) ± jvg''' (ti) satisfies the Cauchy-Riemann equations in ityi > [38], 



,(1) 



and are denoted analytic and anti-analytic (-) functions, respectively. As y ^ 0^ kl^\ti) —>■ g^^^{xi) = 
g{xi) ± jHg{xi), and the limit of an analytic function is an analytic signal. Note that g^^''{ti), i.e. the analytic 
signal evaluated at the complex argument, is still an analytic function for yi > 0, see also [23] [p. 5]. More 
importantly the analytic and anti-analytic signals are represented in polar form by 



(10) 



Any real signal g{xi) G can be decomposed into an analytic and an anti-analytic signal [41]: 



g{xi) 



g^-^Hxi) + g^-Hxi) 



g{xi) + jg^-^Hxi) +g{xi) - 3g^^\xi 



(11) 



= |ag(xi)| cos(27r(/)g(xi)), 

and so a local magnitude ag{xi) and phase <t>g{xi) can for any real- valued function g{xi) be defined from the 
modulus and phase of g^^\xi), whose forms are given by equation (fTOb . We interpret equation (fTTt as the analytic 
and anti-analytic decompositions of g{xi), and this decomposition is important to derive the properties of the 1-D 
analytic wavelet transform, as demonstrated in [56] [p. 426]. These well-known properties of the analytic signal 
need to be suitably extended to the multi-dimensional choice of 'analytic'. 



B. Hype ranaly tic ity 

In higher dimensions to determine equivalents of ag{xi) and cl)g{xi) from an observed real- valued function g{-x.) 
the equivalent/s of function g^^\xi) must be defined and calculated. Based on these definitions the equivalents of 
ag{xi) and (/>g(xi) will then be given some suitable interpretation. For this purpose d-D extensions of the Cauchy- 
Riemann equations are used, and any given extension corresponding to a set of equations wiU be denoted a set of 
Generalized Cauchy-Riemann equations. In 1-D the spatial variable is 1-D, and there is a single associated auxiliary 
variable yi, where an analytic function satisfies the Cauchy-Riemann equations in yi > 0. We consider, in 2-D 
and higher dimensions, augmenting the spatial variable x by a set of auxiliary variables collected in vector-valued 
variable y, of dimension p (some restrictions apply to the choice of y, and its dimensionality), where p will be 
related to the local structure of the image to be represented. We define F^^^ = {y, yj > 0, i = 1, . . . ,p} , and 
consider hyperanalyticity for y G F^*'), and x G M'^. Examples of such spaces include p = \ with yi > and 
X G M*^, as well as x and y of the same dimension, and t = x + jy restricted to the cross-product of d Euclidean 
upper half-planes, usually referred to as a tube see [38][p. 90], denoted Tp'^^ 

We shall be quite careful to distinguish between hyperanalytic functions, and their limits as y ^ 0"*". For future 
reference we make the following set of definition: 

Definition 2.1: The hyperanalytic function. 
Any vector-valued function in spatial variable x, with associated p dimensional auxiliary variable y, that satisfies 
a given d-D generalisation of the Cauchy-Riemann equations for y G F^*') is denoted a hyperanalytic function. 
In analogue to the 1-D analytic signal, we define the hyperanalytic signal. 

Definition 2.2: The hyperanalytic signal. 
Any vector-valued function denoted g^^ ^\^): that can be written as the limit as y — > 0+ of a hyperanalytic 
function kg^ ^^(x, y), with p-dimensional auxiliary variable y is denoted a hyperanalytic signal. 
^(H ^)(x) will have p number of +'s to denote the dimensionality of the auxiliary variable y. In 1-D we found that 
evaluating g^^\xi) at xi taking a complex value ti, the analytic signal still corresponded to an analytic fimction 
in either the upper half plane. This will not in general be the case in higher dimensions and it will not be natural 
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to represent the spatial variable as a complex- valued quantity: furthermore in some instances the dimension of x 
and y will not be the same. For future reference, to treat the rotation of the CWT, we also define a hyperanalytic 
function and signal in a rotated frame of reference. 

Definition 2.3: The 0-hyperanalytic function. 
Any vector-valued function in spatial variable x' = r_9X, with associated p dimensional auxiliary variables y, that 
satisfies a given d-D generalisation of the Cauchy-Riemann equations in variable x' when y is in r(p)((9) is denoted 
a 6'-hyperanalytic function. 

The 0-hyperanalytic function, once the limit y ^ 0+ is taken, yields a 6'-hyperanalytic signal. Note that the choice 
of y may depend on 9. The suitable choice of hyperanalytic to represent local image phenomena depends on the 
local dimensionality of the analysed image. To be able to make a decision on a suitable local representation, the 
local dimensionality of a narrowband oscillatory signal deserves some further discussion. 

C. 2-D Oscillations & narrowband signals 

In 2-D, assuming oscillations associated with a given period are observed, the local dimensionality of the image 
may be of many different forms. Define fo((/>o) = [/o cos((/)o) /o sin{4>Q)]^ , and note that, if components are only 
present with period l//o, this would imply the signal takes the form: 

r2TT 

5r(x) = / ag(0o) cos(27rx^fo((/)o))#o- (12) 
Jo 

Thus the general representation of the image corresponds to a continuum of directions, and ag(0o) is non-zero 
for many (pQ £ [0, 2it) . Often locally the variation of the image is far more structured than the general form of 
equation (fT2b and admits a much simplified representation. For example, typical patterns such as plane waves would 
correspond to a single direction in the signal, or with a slight generalization of slow spatial modulation i.e.: 

= ag(</>o(x),x)cos(27rx^fo(</>o(x),x)), (13) 

whilst textured features [57] often take the form of an aggregation of plane waves: 

L 

5(x) = ^a3(0?)(x),x)cos(27rx^fo(4'^(x),x)). (14) 
1=1 

Not infrequently we find that L = 2, and if this is the case then with ag((/)Q^''(x), x) = ag((/)Q^''(x), x) and 

(2) (2) T 

(x) = '^(x) + 7r/2 and writing x = x [cos(x) siii(x)] , we obtain 

5(x) = 2ag(</>[|^''(x),x)cos(/oxcos(x + '?5<o^^(x)))cos(/o(x)xsin(x + 0o^'')), (15) 

and the signal is separable. Key local characterising features of the image will correspond to L, the local dimensional- 
ity of the oscillation at spatial point x, /o(x), the local period, and the form of the L pairs |ag((/)Q^(x), x), </)'^') (x)|, 
that is the magnitude and direction of the Ith component. We consider three special forms of local structure: the 
local plane wave, that is L = 1, local separable behaviour, that is L = 2 and </)q (x) = (/)q^^(x) + tt/2, and L > 1 
with no special relationship between 4''^\-k) and (/)g^^(x). In the latter case there is no local separable structure 
present. 

The analysis problem then corresponds to, if confronted with a single real image that naturally fits into the model 
of equation il3l . or for that matter equation il4l . determining ag^(0o(x),x), /o(x), and (I)q\'x.) from the single 
observed image, or equivalently defining the 2-D set of functions playing the role of g^^\xi), from g{x). Once 
this has been achieved, the characterising parameters may be determined from g^'^\x). In analogue with 1-D, a 
natural stalling point for such determination, is a sinusoidal decomposition. 
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D. Quatemionic Fourier Transforms 

The FT decomposes structure into plane waves of exp(27rjif-^x), [18] and does not represent separable structures 
as a single coefficient. A natural tool for analysis of separable oscillations is instead the Quaternionic Fourier 
Transform (QFT) [18]. Instead of using complex numbers, the QFT is defined in terms of quaternionic units [58], 
[59]. An arbitrary quaternionic object takes the form q = qi + 92* + + Qik G EI where G M, z = 1, . . . 4 and 

= = = ijk = —1, whilst ij = —ji = k, ik = —ki = —j, and finally jk = —kj = i. q has conjugate 
q* = qi — q2'i' — qsj — 94^) and q can be written as q = 3ft {g} + I^u{q} , where 3ft{(?} = ^[q + q*] , is the real 
part and {q} = ^[q — q*] , is the pure quaternion part. The QFT is then defined in 2-D by [60]: 

Gqiq) = Tq {(/} (q) = / exp(-27rixigi)5r(x) exp(-27rjx2g2) d^x, (16) 

and inverted by 

5r(x) = {Gq) (x) = / exp(27rij;igi)GQ(q) exp(27rjx2g2) c^^x. (17) 

The QFT can recover separate parity structure in xi and X2 as it separately records four real values at each 
quaternion frequency q, thus forming a natural analysis tool for structure separable over x\ and X2 in terms of 

cos(27r3;i(7i) 003(271x2(72)5 cos(27rxi(7i) sin(27rx2(72)5 sin(27rxigi) cos(27rx2(72)i and sin(27rxi(7i) sin(27rx2(72)i each 
corresponding to a separable function in x\ and X2, that is either even or odd. A local representation of this form 
would require defining more than a single additional component g^^^ (x) at x as both local rates corresponding to 
q\ and q2 would need determination. 

The Fourier transform of a signal (^(x) can also be defined in an arbitrary pure unit quaternion e = iei+je2+fee3, 
where e| = 1, can also be calculated, and is given by: 

/CXD /'OO /"OO /"OO 

/ 5(x)e-2-'=f"^ d^x, 5(x) = / / Ge(f)e2--f"- dH. (18) 

We will refer to the above transformation (fTSl as the Unit Quaternion Fourier Transform (UQFT), but it is also 
known as the Type 3 QFT, see [60][p. 2784]. The UQFT is still interpretable as a sinusoidal decomposition of 
structure, as any pure unit e satisfies the De Moivre relationships [61]. The representation is now in terms of plane 
waves, i.e. cos(27rf^x) and sin(27rf^x), and would require more than a single component g^^\x) as both the 
direction and period of f would require determination from the signal. For signals with simple local dimensionality, 
the QFT and the UQFT form natural representations of the signal. Using the FT allowed us to obtain the analytic 
signal/function representation of c(xi) via the stationary phase approximation and we shaU discuss the hyperanalytic 
signals/functions naturally linked with the QFT and UQFT. 

E. The hypercomplex signal 

To extend the analytic signal to 2-D von Billow and Sommer [17] used partial HTs [27] in both of the spatial 
directions, thus defining the hypercomplex signal. The partial HT of g{-x) in either of x/, 1 = 1, 2, when the second 
argument is fixed, is denoted Tii {g} (x), 1 = 1, 2 and the total Hilbert transform is formed by Ht {g} (x) = 
^2 {'Hi {5}} (x). Ht {g} (x) corresponds to the consecutive operation of two partial HTs [27]. 

Definition 2.4: The hypercomplex signal. 
The hypercomplex signal extension of an image g{l^) is defined by 

5(++)(x) = -HCgi^) = g{^)+ig%\^)+3gf{^) + kg%\^), (19) 

where 

g%\^)=ni{g]{^), g%\^)=n2{g]{^) (x) = {5} (x). 
To confirm that the hypercomplex signal corresponds to a valid hyperanalytic signal we firstly need to define a 
hyperanalytic function and then introduce two auxiliary variables, yi and 7/2- We define the set of functions 

Ug-H{y^,y) =p{xi,yi) *p{x2,y2) * g{^), i'^^]^(x,y) = qi{xi,yi) *p(x2,y2) * g{^) 

(2) (3) 

'"g-M^^y^ =p(2;i,2/i) *9i(2;2,2/2) *fl'(x), t'g;i/(x,y) = gi(xi,yi) * gi(x2,y2) * 5(x), 
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where p{x, y) and qi{x, y) are given by equation 0. Collect the spatial and auxiliary variables in the vector-valued 
variable t = x + j {kyi 7/2) ^- Define the quatemionic (or vector-valued) function: 

A:(+^)(x,y) = 7.,;j^(x,y)+^;j,(x,y)+jz;J^],(x,y)+^^],(x,y). (20) 

It is convenient to express A;^^^(x, y) in terms of its Cayley-Dickson forms [62][p. 1952]. We define Ug^^{-x.,y) = 

%;//(x,y) + jt;J^]^(x,y), i;J*5^(x,y) = v^J]^(x,y) + jz;J^J^(x, y), M^^j:^(x,y) = Ug-ni^^y) + jt'J|J^(x, y) and 

^g;l/(X)y) = ^'g;lf(x,y) + jfi'g;l^(x,y). The functions u|^y^(x,y) and ^^'.^(xjy), with • = i, j, are the simplex 
and perplex respectively, of the Cayley-Dickson forms. Thus it in fact follows that 

kitH\^^y) = 4!H(x,y) +^fH(x,y) = n%{^,y) +3v%{-^,y). (21) 

u^g^jj{t), vl^^fj{t) G C and the pair of functions satisfies the Cauchy-Riemann equations in ti for fixed values of 

t2, in the upper half plane of yi > 0, whilst u^J.'^^{t) S C, v^J.^^{t) G C, and the pair satisfies the Cauchy-Riemann 
equations in t2 for fixed values of ti, in the upper half plane of 2/2 > 0. The hypercomplex system of equations is 
defined as: 

d4!H(t) , dv%(t) _ du%{t) dv%[t) _^ ^^^^ 

dyi dxi ' 9x1 dyi 

du%{t) ^ dv%{t) ^ du%{t) dv%{t) ^ 
dy2 dx2 ' 9x2 9y2 

This set of equations provide a valid 2-D generalisation of the Cauchy-Riemann equations, and thus, as the set 
of functions constituting A;^^^(x, y) satisfy equations Mil for (x, y) taking values in the tube Tr, we may 
denote the latter function hyperanalytic, cf definition (12. U . Clearly as y ^ 0"*", Ug-ui^) — > ^(x:), v'^g^^i^) 

5S^(x), ^ 5S/'^(x), and v'^^J.it) ^ 5S^(x), and thus limy^o+ A:J+^^^(x,y) - g^-'^H^), and so g^++H^) 

corresponds to a hyperanalytic signal (cf definition (12.21 ). 

Furthermore, the hypercomplex signal g^~^^\x), when viewed in the quatemionic frequency domain [17][p. 2849], 
satisfies G^^^(q) =0 Vq ^ this providing a 2-D version of equation (0. 

To be able to consider separable structures in an arbitrary direction we extend the definition of a hypercomplex 
signal to a 0-hypercomplex signal. 

Definition 2.5: The 0-hypercomplex signal. 
The 0-hypercomplex signal of a real-valued image g{-x.) is defined for any 9 G [0, 2it) , as 

5^++^ (x) = nCggi^) = <7(r_ex) + ig^,^ (r.g^) + jg^ {r^gx) + kg^ (r_ex) 
= Rgg{^) + iRgHigi^) + jRen2g{^) + kRgnin2g{^) 

= 5e(x)+i5j;i(x)+J5j;i(x) + %i'i(x), (23) 

thus defining the functions 5g*]^(x), for s = 1, 2, 3. The 6'-hypercomplex signal is suitable for analysing structure 
that is separable, or aligned with variations, in a rotated frame of reference. To be able to compactly note the 
QFT of (7g^^^(x) define the reflection operator J'g. The reflection of g((x) in the line that makes an angle 6 with 
the xi axis is given by Jggi'x.) = g{Jgii.), using matrix Jg = [{cos(20), sin(20)} , {sin(20), — 003(26*)}] , where 
6 G [0, 2-k) . Then: 

F[g^^^'^]{i) = [l + sgn([r_ef]i)][l-fcsgn([r_ef]2)]G(r_ef) (24) 
■^Q{5r^^}(q) = ^^[l + sgn([r_,q]i)][l+sgn([r_,q]2)]G(r-,q) 



+ ^4^[l+sgn(Mi)] [l+sgn([req]2)]G(r_0j,/2q) (25) 

^Formally ti = xi + iyi, but any number system of the form of a real number plus a scaled pure unit quaternion is isomorphic to C, 
and so we will relax this formality and not distinguish C = {a:: + jy} from {x + iy} . 
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The result follows by Appendix II-AI combined with the fact that the QFT of a real-valued signal ^(x) can be 
calculated from its Fourier transform (see Hahn & Snopek [26] as well as Pei, Ding and Chang [60]). Thus the 
FT of gg^^\-x.) is supported on [r_ef]2 > 0, whilst as noted in the appendix, the QFT if implemented in a rotated 
frame is supported wholly on one quadrant. Most importantly, as the QFT can be inverted into the spatial domain 
signal, if a signal admits a QFT representation of equation iA-ll then it is the QFT of a 0-hypercomplex signal 
identified by the form of G(r_0q). This result can be used to interpret the properties of the quaternionic wavelet 
coefficients. 

We shall define analysis methods suitable for signals observed in an alternate frame of reference, as the local 
orientation of a signal is a very important local characteristic. Define the signal g_5)(x) via 

5-e(x) = g{r0x), g{x) = g_g{r_gx), 9 G [0,27r), (26) 
sothat5i^J(x) = 5_0(r_0x)+i5i^^.j:^(r_ex)+j5i^^.j^(r_0x)+%i^].^(r_0x), org^^g^J{x.) = g{-K)+ig^^l.^{r_eyL) 
+ jg-l-ni'^-e^) + >^9-l,H(^-e^)- course, g_e{r_gyL) = ReR_eg{x.) = c/(x), but g^fl^e-n = d-l-A'^-e^) = 
RgHsR-egi^) / 5h^(x), s = I, 2, etc. Hence g_Q{r_0:si.) corresponds to the observed signal, whilst 'gilg.j^ 
corresponds to the HTs in a rotated frame of reference observed in our frame of reference. 

Definition 2.6: ^-hypercomplex: hypercomplex vs anti-hypercomplex decomposition. 
The decomposition of real- valued g{-x.) into four quaternionic components ^^^•'(x), ^^^^^(x), ^*^.g^(x), and 

5(x) = i(5i++)(x)+gi-J)(x)+5VW+^iM^ (27) 

where we define ^^^^(x) = 5(x) + i]filj^{r_gx.) + jg':^ljj{r^gx.) + kg''^].^{r_e^), g^^Sg^J {^) = g{^) - 

i9%Hi-r-e^)+39-l-Hi^~e^)-kg%Hi^-e^),l^^^^^ = 9{:!^)+ilfLl.Hi^-e^)-j&,H(.'^-e^)-'^9%H(.^-e^) 
and ^^g.g^(x) = (7(x) — i]j^^^Q.^{r_0-x) — jg^^l.jj{r_QX.) + fe(7^g.^(r_gx). will be denoted as the ^-Hypercomplex 

analytic/anti-analytic decomposition of g{'x.). The four components g^^^(x), g^'g^^(x), ^^^^(x), and ^^^^''(x), 
are denoted as the 0-hypercomplex signal, the 6 first anti-hypercomplex signal, the 6 second anti-hypercomplex 
signal and the 9 third anti-hypercomplex signal, respectively. 

In the rotated frame of reference, equation illl provides a 2-D analogue of the analytic/anti-analytic decomposition 
of a 1-D real- valued signal [41]. We may use the 0-hypercomplex signal to decompose any real image g{-x.) into four 
components, where each component in a rotated frame of reference is analytic/anti-analytic in each spatial variable 
separately. We characterize local structure in terms of magnitude and phase of ]fl^g^J (x) , interpretable in the rotated 

frame of reference. Any of g^J^g^J {x.) if calculated at all 9 = [0,7r/2] , will enable us to parameterise separable 
oscillations. The 6'-hypercomplex signal is introduced in this article for the purpose of recognizing separable or 
univariate local structure that is not aligned with the axis of observation; the act of rotations is integral to defining 
the CWT in subsequent sections, and so its introduction is a necessity. 



F. Properties of the hypercomplex signal 

A hypercomplex signal can be decomposed in terms of a modulus and a set of phase/s, using quaternionic polar 
coordinates, see discussion by [63]. The modulus is still a real number and can be interpreted as local energy content 
whereas the phase/s are given by a unit energy quaternion and contain additional information about the structure 
of the image. This, extended to the 0-hypercomplex signal, corresponds to a polar representation of variation in 
orientation 9. 

Definition 2. 7: The polar representation of the ^-hypercomplex signal. 
The 0-hypercomplex image will be represented in terms of its amplitude, and phases. Define the unit quaternion 



by.r(x) = 'zj^^( 

corresponds to: 



X 



+ iqf\y.)+3qf\:ii)+kqf\yi) = e2™«(x)e27rfc7«(x)g27rj/3,(x)^ ^^^^ ^^^^^ representation 



,(++), 



X 



(+- 
9e 



X 



g27riae (x) ^irk^ie (x) ^2n jf3e (x) 



(28) 
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Von Biilow and Sommer considered only the hypercomplex signal (i.e. 6 = 0) but the exact algebraic forms of the 
angles for both ag{x.) and f3e{'^) in terms of 51^ ^^(x), 1 = 1, 2, 3, can be found by substituting in the form for 
the quaternion components into the formulae given by [17] [p. 2849]: 



1 

— tan 
47r 



1 

— tan 
47r 



1 

47r ■ 



■ sm 




,(3)/ 



x)gf (x)+g^^^(x)5^^^(x 



,«^x)+,f)^x) 



(x)g;^^(x)+g;^^(x)gr(x 



.(2)2 



(x)+gr(x 



(1)2, 



X + 



.{2)2, 



X 



.(3)2 



(x)+gf'(x) 




X 



X 



(29) 



More than one phase function are used, as variations in either axes are given by ag{'x.) and /36i(x) respectively, 
whilst 70 (x) is considered by von Biilow and Sommer as a measure of 'degree of separability'. If the signal 
analysed takes the form g((x) = gi{xi)g2{x2), i.e. the signal is separable, then they note that [17][p. 2849]: 



ao(x) = tan 



1 f m{9}(K) 

9(x) 



/3o(x) = tan 



1 f H2{9}{^) 
9(x) 



, and 7o(x) = 0. Even if (^(x) is not separable the polar 



representation of equation d28b may be used to represent an arbitrary real- valued image ^(x). 

Definition 2.8: The 6'-polar representations of g(x). 
The 0-polar representation of a real signal g{x.) is given by: 



5(x) 



5R- 



g27ria _ e ,e (x) g27rfc7_ e ,e (x) ^-KjP^e.e (x) j_ 



51'^ ^' (xj (^cos(27rQ!_e,6i(x)) cos (2717^0,61 (x)) cos(27r/3_e,e(x)) 
+ sin(27r5_e),6i(x))sin(27r7_e),6i(x))sin(27r/3_6i,e(x))^ . 
A separable structure in orientation 9 will have 7^ (x) = and the signal is represented by: 

5(x) = 



(30) 



(31) 



cos(27ra_0^5)(x)) cos(27r/3_g^0(x)). 
Note that S-e.elx) etc are determined by equation (l29t but where now the ^-hypercomplex extension of ^-^(x) 
is used as the basis of the representation. To see the utility of such a representation consider analysis of a signal 
separable in the rotated frame of reference, i.e. 



9W 



(32) 



Under the assumption that g{-x.) satisfies equation i32l . we determine that the 0-hypercomplex extension of the 
signal is given (in terms of x' = r_ex ) by: 



~(+4 

9-e4 



5_e(x') + iHi {g-e} (x') + {g-e} (x) + ^^2^1 {^-4 (x). 



(33) 



TCi {g-e} (x) = Ti {gi^-e} ixi)g2,-eix2), 

'H2{g-e}{^) = T-i {g2-e} {x2)gi-0{xi) and ^2^1 {^-e} (x) = 'H{g2-e} {x2)'H{gi-e} (xi)- Hence, in our 
extended definition, if (^(x) G M is separable in some rotated frame of reference x' = r.gx, then 7_e,6i(x) = for 
some e G [0,7r/2). 

If a single component is present in g{-K) at x and the aim of the analysis was to determine the orientation of the 
separable structure in the image at this location, denoted 9', then we calculate g^^g^J {-x.) at all angles 9 € [0,7r/2), 
and determine the value of 9' by 7_£)' g/(x) = 0. Once this has been done the signal is well described by its local 
magnitude, i.e. 



whilst local variation is best described by the two separable phase functions 

5_e,e(x) = tan"^ {n{gi-e} {x'i)/gi-e{x'i)) , /3_e,e{^) = ^an""^ {52,-0} (4)/52,-0(4)) 
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The local frequency is found in each of the two axes as -^a-e^g{x.) and gp-/3_6i,0(x). 

Proposition 1: Orthogonality of the signal and the HTs. 
The original signal g{x) and g^^\x.), s = 1,2,3, are all mutually orthogonal. 

Proof: This is a trivial consequence of the orthogonality of a 1-D HT. Furthermore we note that as g_0{x.) is 
orthogonal to ^^^^.^^^(x), s = 1,2,3, we find that g{x.) is orthogonal to 'gfl.^{r^gx.), s = 1,2,3. 

Proposition 2: Norm of the hypercomplex components. 
The g^^\ s = 1,2, 3, have the same norm as g{'x.). 

Proof: This is a direct result following from the norm of the 1-D HT of a real-valued function [23] [p. 4]. 
Hence in summary: {g^'^^^x.)^ defines a quaternionic object where each real-valued function multiplying the 
quaternions i, j and k can be considered as objects with the same spatial/spatial frequency support as ^(x), as 
the HT constructs an object with the same time-frequency support as the original signal [64]. Note that all of 
|5ff0 6»(-^)} s = 1, 2, T are mutually orthogonal. ,f^u |(7i^^^(x)| constructs a quaternionic object with the 
same spatial and spatial frequency support as ^(x), where the relative magnitudes of the components specifies the 
structure of g(x) in a rotated frame. 



UB.;gi^,y) = ipi-,y)**9{-))i^), v'-2(x,y) = {qi**g){x,y), (35) 



G. The monogenic signal 

The hypercomplex signal, once extended to the 0-hypercomplex signals recognizes locally separable structure, 
and can if the signal is calculated for all 9 E [0, 27r) , determine locally unidirectional structure at a given orientation. 
However, this does require the calculation of the hypercomplex extension at all 6. A hyperanalytic extension for 
locally unidirectional structure, or plane waves, at any given spatial point can also be made. Felsberg & Sommer, 
[25], introduced a hyperanalytic extension in image processing suitable for this purpose. See also an early work by 
Dixon [39], where monogenic functions were denoted 'pure Hamiltonian' functions and the discussion by Duffin 
in [65]. This approach starts from the d-D Poisson kernels [38], a d-D version of equation Q, defined for any 
dimension d> I : 

y xi r(^) 

p{x,y) = Cd qi{x,y) = Cd j+r Cd = — I = I, . . . ,d. (34) 

[|x|2 + 2 [|x|2 + y2] 2 TT 2 

An extra set of d + 1 functions 

are defined from g{x). For d = 2, ^UR-g{x,y),v^^^g{x,y),v^^^g{x,y)^ satisfy an alternative set of generalized 
Cauchy-Riemann equations, with only a single auxiliary variable y. These generalized Cauchy-Riemann equations 
are called the Riesz system, see [38][p. 234]. 

Definition 2.9: The monogenic function. 
Any solution of the Riesz system in the upper half-space (y > 0) is denoted a monogenic function [66] [p. 35]. 
Note that if k^^\x,y) is a monogenic function then A;(+)*(x, y) is a solution of the Riesz system in the lower 
half-space {y < 0), and will be denoted an anti-monogenic function, with notation k^^\x,y). 
For d = 2, if k^^.]^{x,y) = UR-g{x,y) + if^.^^(x, y) + jt;^?^(x, y) then it is a monogenic function [38][p. 235]. In 
the limit of y ^ 0+ [38]: 

lim UR.g{x,y) = g{x), lim v'^^ {x,y) = g{x) * ...*ri{x) = 7^^5^(x) = g'-^ix), (36) 
J/— >o+ y^o+ 

where the quantities n(x) and Ri{f) are given by: 

ri=CnY^, Rl{i) = -j^, Ri,Q{q) = -i^, i22,Q(f) = -j- ^ = 1,2. (37) 
\^\ J Q Q 

{r/(x)} are the Riesz kernels for d> 2, and the FT and QFTs of the Riesz kernels are given by [27], [38]. 

Definition 2.10: The Riesz transform. 
The Riesz transform of a signal g{x) € L^(M^) is defined by operator TZ 

n{g} (x) = in{g} (x) + j7^2 {g} (x), TZi {g} (x) = ^^^(x) = (n * *g) (x), / = 1, 2. (38) 
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Felsberg & Sommer [25] recently introduced the monogenic signal into image processing by combining g{-x.) with 
n{g} (x): 

Definition 2.11: The monogenic signal. 
The monogenic, and anti-monogenic signal of real signal ^(x), are defined by applying operator to signal 
5(x): 



<7(±)(x) = >^±M(x)=5(x)±7^M(x) = <7(x)± ig'^>{^) +jg'^>{^) 



(39) 



Note that TZ, and g^^ (x) are defined by equation d38l) . 

Letting y ^ 0^ the monogenic signal is retrieved from the monogenic function, and so we find that (/(+Hx) IS 
a hyperanalytic signal, c/ definition (I2.2I) . See also the careful discussion in [66] [p. 32-37] on monogenic signals 
and the generalized Cauchy-Riemann equations. 

Only one auxiliary variable is introduced when constructing the Riesz components, linked with the norm of x, 
and the monogenic and anti-monogenic signals play the roles in 2-D of the analytic and anti-analytic signals in 
1-D. To accommodate that the observed image is not observed in the direction of the local variation, we introduce 
here the concept of a ^-monogenic and anti-monogenic signal. 

Definition 2.12: The ^-monogenic signal. 
The ^-monogenic, and 0-anti-monogenic signal of real signal g(x) are for any value of a rotation 6 G [0, 27r) given 
by: 

5^^) (x) = ReM^ {g} (x) = {g} (x) = g{r_e^) ± [iReUig^^) + 3Ren2g{^)] (40) 



g(r_ex)± ig''^\r_e^) + jg^^\r_g^) = ge{^) ± ig^!^(x) + jf5r^:'jj(x) 



thus defining the components <^g^]j(x) for s = 1, 2. 
Note that 



J' {g'i\r^e^)} = -j cos((/. - 0)G(r_,f), {^^^(r^.x)} = -j sin(</. - 0)G(r_,f). (41) 

Theorem 1: ^-monogenic functions and their monogenicity. 
The ^-monogenic and the ^-anti-monogenic signals g^^\-x) are both the limits of sets of functions satisfying the 
Riesz system of equations in a rotated frame of reference. 



Proof: See section II-BI Thus a 6'-monogenic signal can be considered to be the limit of a monogenic function 
in the rotated frame of reference. For consistency with our notation so far, if 6* = 0, the signal is denoted g^^\-) 
without subscript. Furthermore we may note that 

G^^) = G(r_ef)±(-fccos((/.-0)G(r_0f) + sin(0-6')G(r„ef)) 

= [1 ± sin(0 -e)^k cos{(p - 6)] G(r_ef ). (42) 

If the Fourier transform of a function can be written in the form of equation (l42t it can be directly deduced that the 
function corresponds to the ^-monogenic image of ^(x). Thus just like we intend to use equation (E^ to recognize 
a quaternionic object as corresponding to a 0-hypercomplex object, equation d42t can be used to determine if a 
quaternionic object is ^-monogenic, and what real- valued object plays the role of ^(x). 

A real- valued image can be decomposed into a ^-monogenic and 6'-anti-monogenic component. Using 50 (x) 
defined in equation (l26t we represent g{-x.) by: 

= 9-eir-e^) ± (^ig^jiir^e^) + jg^^lj^ir^gx)^ 

= g{x)±(i^}l.^{r_e^)+3g^%,^{r_e^)y (43) 
TZiR^gg{x) = cos(6')i^_e7^l5((x) + sin(6')i^_e7^2S'(x), and as noted by [38][p. 241]: 

5i^j.^(r_ex) = ReUiR^egi^) = cos(e)7^l5(x) + sin(e)7^25(x) ^ nig{x) 

g^Sni^-e^) = ReH^R-egi^) = - sin(e)7^l5(x) + cos(0)7e25(x) / 7^25(x). (44) 
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The rotation operator, Rg, and the Riesz transform operators, TZi, 

7(±), 



1 = 1, 2, do not commute, and so 



(x) / 



5^^'lxj. 

Definition 2.13: ^-monogenic: monogenic/anti-monogenic decomposition. 
The decomposition of real-valued g{x.) into two quaternionic components g^^Jgix.) and 'g^_} J'x) : 



9W 



(+) 



(45) 



will be denoted the 0-monogenic/anti-monogenic decomposition of ^(x), where the two components are given by: 



^(x) lb ^i5r^^g.j:j(r_0x) + jg^_^0.ji{r-e^)j • This is a natural decomposition of a locally 1-D signal more 
naturally observed in a rotated frame of reference. Equation ( R3t forms a 2-D extension of the decomposition of 
a real- valued object into a quaternionic form in analogue with the analytic/anti-analytic decomposition of a 1-D 
signal, see [41]. 



(2) 



H. Properties of the monogenic signal 

Again, it is natural to represent a hyperanalytic function in terms of amplitude and local structure, and so [47], 
[66], [67] represent the monogenic signal in terms of an amplitude and a phase. A single phase only is used to 
describe variation as the monogenic representation corresponds to a local plane wave structure. The phase represents 
the local period in the direction of variation determined from the signal. 

Definition 2. 14: The polar representation of the ^-monogenic signal. 
The ^-monogenic and anti-monogenic images will be represented in terms of the amplitude, phase and direction 
by 



(x) = |55+)(x)| e±2-.,(x)<^.(x)^ 
where the amplitude, phase and direction are given by 



[5e(x)]2 _^ [5fJ^jj(x)]2 + [5^f^(x)]2, , e^,(x) = [i cosue{x.)+j sinz^e(x)] , 



.(2) 



(46) 



(47) 



tan 



9iU^) 



1 

X) = — tan 
^ 27r 



+ bS(x)]2 



The sgn(-) added to the definition of (/>e(x) to appropriately determine the range of 
and 0-anti-monogenic images of a real image g{x.) have the same modulus and phase. g^g^^*{x.) 

Definition 2.15: The ^-polar representation of g(x). 
The ^-polar representation of a real signal g(x) is given by: 



5e(x) 

x) . Note that the ^-monogenic 



,(-), 



X 



5(x) 

Thus locally a signal is described by its magnitude 



~{+) 

9-0.6 



:x) 



and the local structure of the signal is determined 



by 5i(x), from which an instantaneous frequency is determined by V0_e,6»(x) • From equation (I44t it clearly 



follows that: 



9%{^) 



\9 



(±)( 



X 



but ei>_,^,(x) 



/ ejy(x). The monogenic representation of 



g{x.) is thus equivalent to the ^-monogenic representation, unlike the case of the hypercomplex and 0-hypercomplex 
representation, that are clearly distinct. The intuitive understanding of these relationships is that the monogenic signal 
is constructing a plane wave representation at spatial point x, where the direction of the plane wave is determined 
from the signal. No matter what axes we construct, the local plane wave representation in the amplitude and period 
(phase) remain the same, but the parameterization of the direction changes with the axes. The reason for introducing 
the ^-monogenic signal is that the calculations of the properties of the CWT are simplified. Hence for convenience 



the notation 



g'Sg ei^) (p-e.ei^), are retained, even if both are in some sense redundant. 
Proposition 3: Orthogonality of the signal and the Riesz transform. 
A real image g{-x.) is orthogonal to each component of its Riesz transform TZ. However, the two components of the 
Riesz transform, g^j^\x.), g^^\-K) are in general not orthogonal. 
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Proof: This is a direct generalisation of the n = 1 case where the HT is orthogonal to the original signal. Indeed, 
exploiting the Hermitian symmetry of the FT of a real image G*(f) = G(— f), we have 



{g^^\g) = I dH j ^ G*(f) G(f) = J d'f j ^ G(-f) G(f) = 0. 



(48) 



The function G(— f) G{i)/f is even-even (ee) with respect to /i and /2 whereas /i itself is obviously odd- 

(2) 

even (oe). Since the integration is over the entire E2 plane, it is zero. By the same argument {gp^ ,g) = and 
therefore {IZg^g) = 0. Furthermore we note that is orthogonal to ^|^g(x), and thus g((x) is orthogonal to 

5_g(r_ex), s = 1, 2. However, the two components of the Riesz transform itself are not in general orthogonal: 



{g^^a^R) = I ^G(f)G(-f)d2f ^ 0. (49) 



The natural way to view the monogenic image is as consisting of two orthogonal components, (^(x) and Tlg{x.). 

Proposition 4: Orthogonality of the components of the Riesz transform of an isotropic function. 
If the signal g{x.) is radially symmetric and real, then the integrand in ( I49t is 

(5g\52^)=0. (50) 

Proof: G(f) = G{f), and radially symmetric as well as real, whilst /1/2 is odd with respect to either /i or /2. 
Thus {g^^\g^^^) = 0. Note that for an isotropic signal g((x) =^_5i(x) and so g{x.) is orthogonal to ^^g^(x). 

Proposition 5: Norm of the Riesz components. 
The norm of the Riesz components is given by (as noted by Felsberg & Sommer [25] [p. 3140]): 

ia^R^a^R) = I ^G'(-f)G(f)d2f, s = l,2, {g^i\g^i^) + {g^n^,g^n^) = {g,g). (51) 
If furthermore 17 (x) is radially symmetric then 

{g^R\g^n^) = l{g,g), s = 1,2. (52) 

These relationships were proved and used in [44] in the special case of the isotropic multiple Morse wavelets. In 
summary 7^5 (x) is orthogonal to ^(x), has the same frequency support, as well as the same total norm as (^(x) 
and exhibits similar spatial decay as (/(x) (or at least as similar as 7{g{xi) does to g{xi).) Hence as Hg{xi) is 
considered as having the same time-frequency structure presence as g{xi) (see for example Boashash [64]) and we 
consider TZg{'x.) as having the same spatial and spatial frequency support as ^(x). 



/. The monogenic signal and the UQFT 

The monogenic signal determines the directional preference of a given real signal at a given point in space, x. 
Let us discuss this aspect in some more detail. As the monogenic representation is suitable for a plane wave field 
51 (x), consider the representation given to a purely unidirectional signal (7c/ (x). By assuming that 5(7 (x) purely 
experiences variation in a given direction, implies that it can be represented as a superposition of plane waves with 
direction v. Then for some given constant n = [cos(z^) sin(i/)] : 

/•oo 

gui^)= Gc/(/)cos(27r/n^x)(i/. (53) 
Jo 

Define the unit quaternion e,^ = cos{v)i + sm{v)j. We calculate the monogenic signal of gu{^), and determine 
that: 



roo 

g'l^'i^) = / Gt;(/)exp(2^e,/n^x)(i/, 
Jo 



roo roo 

= / Gu{f)cos{27Tfn^^)df + [icosiu)+jsm{iy)] sin(2^/n^x)d/. (54) 

Jo Jo 

This result follows from the monogenic extension of a sinusoid given in [25], and the global orientation of the 

globally directional signal may be determined from the monogenic signal by calculating the ratio of the two Riesz 
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transforms. Let us now demonstrate how gfj (x) is related to the UQFT in dy. It may be shown, again using the 
monogenic extension of a sinusoid given in [25], that g|j*2 (f) = Gjj{f), / > 0, so that the spectral support of 
the UQFT of the signal in is one-sided, and limited to the positive frequencies. Thus an equivalent to equation 
(0 can be derived that corresponds to an alternative to GQ^^^(q) = 0. The result is to be expected, as the plane 
waves are really a 1-D feature embedded in a 2-D space. If the signal corresponds to unidirectional variation only, 
the monogenic signal is extracting the orientation, and lives in a half-plane, whilst the anti-monogenic signal lives 
in the second half-plane. Calculating the partial HT in the correct direction will yield an equivalent representation 
to using the monogenic signal, but the advantage of using the latter is that the direction is retrieved directly from 
the observed image without having to perform HTs at all values of 9. The assumption of global unidirectionahty is 
very strict, and constructing a local description of the signal will be essential for the performance of the analysis 
method. 

III. Wavelet analysis 

Wavelets provide a method of constructing a local decomposition of an observed image, that for many classes of 
images constitute of decomposition coefficients mainly of negligible magnitude. Furthermore, due to the definition 
of the wavelet coefficients, it is also frequently easier to interpret the non-negligible coefficients, rather than trying 
to disentangle the full behaviour of an image in the original spatial domain. The interpretation of the coefficients 
relies on the choice of mother wavelet. To be able to satisfy the reconstruction formula of equation Q for any 
arbitrary function (^(x) G L^(M^), the mother wavelet function, denoted V'(x), must satisfy the two conditions of 

< < oo, where = {lirf j ^^^^^ d^f, j |V'(x)|^ d^x = 1. (55) 

Thus the mother wavelet is constrained to be both a spatially local function, and to be mainly supported over 
some range of frequencies not including the origin, i.e. the function is oscillatory and zero-mean. Without loss 
of generality, center the mother wavelet function in space to x = 0, and assume that ^'(f) is maximum at f = 
fo = /o [cos((?!)o) sin((/)o)]^ • Of course the possible structure V(x) may exhibit whilst still satisfying equation 
d55b is very variable, and even more so in 2-D compared to 1-D. For example, the wavelets developed in [20], 
[21], [43], [68] are complex and describe directional plane waves, as do the wavelets of [1], whilst a great deal of 
applications that rely on perfect reconstruction use discrete filter bank separable wavelets [4], and a third possible 
option corresponds to using radial functions, as described in Antoine et al. [13], Metikas & Olhede [44] and also 
by Farge [69]. Visually, in a contour plot description of the wavelets, highly directional plane wavelets have the 
appearance of local maxima corresponding to evenly spaced hnes with a given normal, separable wavelets have 
local maxima and minima in a checkerboard patter, and isotropic wavelets correspond to local maxima appearing 
in repeating ring-shapes with approximately an even spacing between the rings. The 'optimal' choice of wavelet 
naturally depends on the application in question. 

This article aims to provide classes of mother wavelet functions so that the wavelet coefficients calculated are 
hyperanalytic signals in the index b. Additionally energy in the observed image should be separated to the correct 
decomposition coefficient, and this corresponds to using wavelets that are essentially supported over a limited 
region of space and spatial frequency. To satisfy the latter restriction wavelets that are well-localised [13], [44], 
are often used. Local structure will be represented in terms of interrelationships between coefficients in a vector- 
valued representation at each local decomposition coefficient index point. To construct a four-vector (quatemionic) 
description of local structure, we use wavelets of the form: 

V'(x) = V'^''^(x) + #(*)(x) + jV^^'^(x) + k^^^\-iL), (56) 

where each -i/'^'*'' £ L^(M^) for s = r, i, j, k. Thus for any x G it follows that '(/'(x) G H. 

In general we construct V'(x) from il)^^\yL) and will denote the construction operator by V(x) = WV'^''^(x). 
If the observed image can locally be described as a plane wave, as is the case in equation ( I13t . then the local 
monogenic description of the signal is appropriate, and the CWT will be constructed to correspond to a local UQFT 
(c/ equation dlBll.) If on the other hand, local separable variations in an appropriate frame of reference are present 
like the signal modelled by equation M5\ . then a local version of the QFT, cf equation (I16t . should be used for 
analysis. 
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Naturally the whole family of wavelet coefficients constructed from the family of wavelet functions must enjoy 
suitable properties. 

Definition 3.1: Members of the quaternionic wavelet family. 
An arbitrary member of the quaternionic wavelet family is defined by 

V'c(x) = V(x) = ^^'^(x) + iV-^^^W + jV^^'^(x) + kijf\x). (57) 
The construction corresponds to individually translating, scaling and rotating the three components of the quater- 
nionic function. To define any member of the quaternionic wavelet family, that is ip^{'x.), each of the real wavelet 
functions that combine to make up the quaternionic wavelet is scaled, rotated and shifted so that ip^i'x.) = 
W tp^'^\x.). Note that in general ?/'^(x) ^ W ^^^\x), and some care must be taken when noting the 
properties of the members of the wavelet family. Also the properties of the CWT coefficients must be considered 
with some care, as in general neither is wy\;^{r) g) = Ww^(r) g), where W is implemented in b. It is very 
important that the wavelet coefficients are interpretable, as they may be considered as a local projection of the 
image, and they, rather than any other given quantity, will be the basis for analysis of a given image. We shall 
discuss the properties enjoyed by the wavelets as well as the CWT coefficients in detail. 

A. Hypercomplexing wavelets 

The QFT decomposes an image into sinusoids, where the parity of variations in both xi and X2 can be 
considered separately via the i, j and k components of the QFT. From equation (fT6t we can note that the 
construction of the QFT is highly non-trivial, as ^(x) is left-multiplied by exp(— 27rixigi) and right-multiplied 
by exp(— 27rjX2(72)i so if a quaternionic g{-x) is decomposed in the QFT coefficients, then the transform can- 
not be described in terms of an inner product of (^(x) with another function. If g(ii) G M then Gq{(\} = 

f^oc di^) exp(-27ri2;igi) exp(-27rjfa;2g2) 
d^x, and the QFT may thus be expressed as an inner product between [exp(— 27rixigi) exp(— 27rjX2(72)]* and ^(x). 
Note that 

(^^-2niq^x^^-2njq,x,y ^ (x) + (x) + jl?(^°)(x) - fcl9(°°)(x), (58) 

thus defining the functions {t?*^'')(x)} , s = ee, oe, eo, oo. = cos (27r (71X1) cos(27rg2a^2) is an even function 

in both xi and X2, = {i?^^'^)} (x), ??M(x) = H2 {i?^^^)} (x), whilst i9(°°)(x) =ni{n2 {'d^"''^}} (x). 

To duplicate the structure of |i?(' )(x)} , constructing a local QFT, a real-valued mother wavelet function ^lJ^^^\^x) 
is first chosen to be even in both xi and X2, but not necessarily separable. We then define the hypercomplexing 
wavelet, by adding three extra functions to this function, where each extra function multiplies a unit quaternion, 
and where the functions are odd instead of even in either (or possibly both) of the two spatial directions. 

Definition 3.2: Hypercomplexing wavelets. 
The hypercomplexing wavelet function ^/;++ ( • ) is formed from a 2-D real valued wavelet function ip^'^^^x.) G (M) , 
satisfying the admissibility condition of dSSl and additionally the condition of: 

V,{-)(xi,X2) = V^'=^)(±Xi,±X2) = V^'=^)(X2,X1), (59) 

i.e. the real wavelet is even in both xi and X2 and symmetric across the indices. For such ^'^'^'^^(x) let: ijj^"'^\x.) = 
Hi {i;^^^^} (x), '0^^°^(x) = W2 {^(^')} (x), and V^^^'^x) = W2 {Wi {V'^'^^}} (x), and form the hypercomplexing 
mother wavelet via 

V'+^Cx) = T^CV'^'^^ (x) = V^^'^ (x) + (x) + jV'^'") (x) - (x) . (60) 

The hypercomplexing wavelet is defined by a set of three functions constructed from the partial HTs of the even/even 
function in the two spatial directions, but note that HC / HC, (where HC is given by equation (I19t.) The Fourier 
transforms of the functions forming ■tp^~^{x.) are given by: 

^(-)(f) = (-j)sgn(/i)M'(-)(f), M/M(f) = (-j)sgn(/2)vI/(-)(f), 

= -sgn(/i)sgn(/2)^(^'=)(f). (61) 

^(ee)(f) £ as V^''^^(x) is even in both xi and X2, and thus ^'^^''(q) = ^ (q). 

M/++*(f) = (l-fcsgn(/i)+sgn(/2) + fcsgn(/i)sgn(/2))*^(^'=)(f) 

= (l + fcsgn(/i)+sgn(/2)-fcsgn(/i)sgn(/2))^'(^^)(f). (62) 
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The Fourier transform of the convolution of t(;~^'^*{-x.) with g{-x.) G (M?) , taking values in M is given by 

T {g * *V^++*} (f) = G(f) (1 + fcsgn (/i) + sgn (/s) - fcsgn (/i) sgn (/s)) ^(^^^(f). 

Proposition 6: Admissibility of the hypercomplexing wavelets. 
The quaternionic function '0+"^(x), as well as any of its four component functions, all satisfy the admissibility 
condition. 

Proof: As ^lJ^'^^\x) is an admissible wavelet function, it follows all its partial HTs are admissible, i.e. square 
integrable and satisfying the admissibility condition given by equation dSSl . 

Theorem 2: Lack of hypercomplexity of tp'^'^ (x) . 
The function tp^^ (x) defined by equation (I60t is not hypercomplex. 

Proof: By direct calculation we note that the QFT of ip'^^ix.) takes the form 

/ if qi<0 or 92 < 0. (63) 

Thus the QFT of V'^^(x) is not wholly supported on positive quaternionic frequencies only, and ^/;++(x) is not 
a hypercomplex signal. The hypercomplexing wavelets are therefore denoted ^p^~^{■), rather than ip^^'^^-), but 
the wavelet coefficients they construct are w^^^\^;g), as they are hypercomplex signals in b. Note that Chan et 
al. [2] define discrete quaternionic wavelet filters by combining two real-valued discrete wavelet functions that are 
approximately HT pairs, in a discrete version of equation (I60t . The objective of the quaternion-valued analysis is to 
form representations of signals that are hyperanalytic signals rather than using hyperanalytic signals in constructing 
the decomposition, and as we shall see, the hyperanalytizing wavelets are useful decomposition filters. 

Proposition 7: FT and QFT of the members of the hypercomplexing wavelet family. 
Any member of the hypercomplexing wavelet family defined by substituting in the form of the mother wavelet into 
equation (1571 has a FT and a QFT given by 

vl/++(f) = \a\ [1 + sgn ([r_ef]2) - fcsgn ([r_,f],) [1 - sgn ([r^^fls)]] ^^''\ar_e^)e-''^i"^^'° 



\a\ e-2-^«^^^ iy^*(-)(ar_eq) (1 + sgn ([r_eq]i) + sgn ([r^^qla) 

1 + 

-sgn([r_0q]i)sgn([r_eq]2)) + ^'('^^)(ar_e J^/2q) (1 + sgn ([r^q] J 

+sgn ([req]2) - sgn ([r,q],) sgn ([r.q]^)) e~^^i''^''^ (64) 



respectively. 



Proof: The Fourier transform follows by direct calculation, whilst substituting in the form of the mother wavelet 
into equation ( 1631 yields the expression for the QFT. Given the mother hypercomplexing wavelet was not a 
hypercomplex signal it is not surprising that neither will an arbitrary member of the family be a hypercomplex 
signal, as we may note from equation (l64l) . We will build hypercomplexing wavelets from either ijj^^'^) (x) = ip^^^'^ (x) , 
the separable mother wavelet corresponding to a tensor product of two 1-D functions, or the isotropic real mother 
wavelet ^'^'^^^(x) = ^p^^^\x), that is a function of x = y^^^i~+~2?2- The reasoning behind these choices, is that we 



replace '&^^^\^) by functions that are even in both spatial directions, and that are chosen to satisfy equation ( I59t 
Either of ip'f^{^) and ip^^{-x) will enable us to decompose the observed image into local 0-hypercomplex signals, 
for any given values of a and 9. 

1 ) Separable hypercomplexing wavelets : The separable hypercomplexing wavelets are constructed from a single 
1-D even function, denoted ^[^(xi). ^'^^(x) is constructed to filter the signal to contributions corresponding to 
the same period in both spatial directions - by augmenting the doubly even function with the three extra wavelet 
functions we construct directional wavelets that are well localised in both direction and scale as well as spatial 
position. Using ipg^{-K) we separate components according to their directional localisation, and the CWT needs 
to be calculate across the fuU range of G [0, 7r/2) . First define the doubly even mother wavelet that will serve 
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Fig. 1. The real part ((a)), first, second and third imaginary parts ((b), (c) and (d)) of the separable hypercomplexing Morse wavelets, with 
ni = n2 = 0, l3 — 9 and 7 = 4. 



as the real component of the quaternionic function 'ip^^{x.), by t/j^'^'^^x.) = ip[^^{xi)ijj[^^{x2)- We define the 1-D 
function ip[°^{x) = H jT/^i*^} (x). We may then observe that using equation (1611) in this instance we find: 

Assume /o = arg^^Q max^^'^(/). The wavelet functions may be considered a filtering the signal to frequencies 
near /i, /2 ~ ±/o. Note that in 1-D analytic and anti-analytic wavelets are defined by combining Vi'd(^i) ^i^^ 
''Aid(^i) ^ complex-valued wavelet function by: 

= i^'ii^ (^i) ± iV'iS (^i) ^S3(/i) = (/i) ± (/i) e M. 

Thus the FT, and the QFT of the members of the separable hypercomplexing wavelet family are given by: 



-27rjf^b 



(65) 



(66) 



(q) 



\a\ e 

1 + k 
+— ^ 



2niqib 



1-k 



respectively. The interpretation of ^'^^(f) is that the hypercomplexing wavelets are naturally treated in a rotated 
frame of reference, whilst from the form of the QFT it is apparent that there is no choice of rotation that will make 
the wavelet ^-hypercomplex. However, as the mother failed to be hypercomplex, this is not surprising. 

Using the separable hypercomplexing wavelets each of the two analytic wavelet functions may be given a polar 
representation, i.e. 



e^''^'''^'/'.io(^\ and thus it transpires that V'5^(x) = ipl^'ixi) ^id (^2) 



a^(x) = /3^(x) = 4']^-l^{x2), and 7^(x) = 0. The separable hypercomplexing wavelets are representing 

local structure by finding the local structural information in xi and X2 separately, and representing this in terms of 
a magnitude and two phases. <^^^i£)(2;) is the local frequency content of the 1-D function, Tp^{xi), and in the two 
separate spatial directions the same frequency behaviour is considered. This is obvious from from Figure ^ where 
the wavelets correspond to two separable oscillations local to x = 0, with the same local period in both spatial 
directions. 

2) Isotropic hypercomplexing wavelets : Basing the hypercomplexing wavelet on an isotropic mother wavelet, 
V'j^'^^(x), wavelet coefficients are produced by filtering the signal to retain contributions with the same period. 
Using the extra components potentially separable structures can be characterised that have the same period present 
at a given spatial location. Figure |2 shows isotropic hypercomplexing wavelets, and how the spatial localisation is 
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radial, whilst each added component is odd in xi and X2, respectively. When an isotropic mother wavelet is used 
then we obtain that: 

These wavelets are able to determine any preference for variation in the direction x = tan~^ {x2/xi) = 0, x = "^72 
and X = respectively. This is clear from the example of an isotropic hypercomplexing wavelet given in Figure 

m 

3) Hypercomplex wavelets : Following work on discrete wavelet filters given in [l][p. 138] a real directional 
wavelet based on the 1-D analytic continuous wavelet filters of section HIl- A. 1 1 will be constructed via first defining 
the 7r/4 directional wavelet function: 

V^(,"/')(x) = VS(X1)VS(X2) - VS(X1)VS(X2). (68) 

V'^^^''(x) is local in frequency to f = ib\/2/o [cos(7r/4) sin(7r/4)]^ . Formally, the above wavelet, does not satisfy 
(l59t . and cannot determine existing separable structure as it only considers structure along the diagonal direction 
and not variation perpendicular to this direction. To align the wavelet with the axis of observation we define 

fe(x) = *<-"V,/,.) = «Sg (^1^) (^) - -^SS (^) *S (^) . m 

Thus (x) will be local in frequency to f = ib\/2/o [1 0]^ . The function is extended to a quaternionic object by 
calculating its total and partial HTs: 

nj f I w X / Xi -X2\ ,(e) / Xi +X2\ , (e) / Xi - X2 \ (o) / Xi + X2\ , 

n2 {^d) (x) = 0, Ut {V'd} (x) = 0. 

Formally by direct calculation: ^/;^^^(x) = V'd(x) + iV'D,2(x). ^^^^^ (yC) is a directional quaternionic wavelet 
function, naturally taking the form of a wavelet suggested by [1]. This wavelet is formed by a similar set of operations 
to the hypercomplexing wavelet but we have given it the notation "0^^^ (x) , as formally, it is also a hypercomplex 
signal. As the third and fourth real-valued components of the quaternionic objects are equivalently zero, and it is 
the non-zero fourth component, which prevents the hypercomplexing wavelet from being a hypercomplex signal. 
V'^^''(x) can be treated as a complex-valued wavelet (see footnote but is not equivalent to the complex 
directional cone wavelets of Antoine et al. [21]. The two components of ^/^^^^(x) are plotted in Figure 0] 



B. Monogenic wavelets 

The hypercomplexing wavelets were constructed so they filter the image as a local version of the QFT in 
orientation 6*, and their usage will (as we shall see) produce wavelet coefficients that are 6'-hypercomplex signals. 
The monogenic signal is locally identifying a dominant orientation of variations and constructing, by adding the 
Riesz transforms of the components to the signal, a polar representation of such variation. The monogenic signal 
of a globally unidirectional signal could be represented as an Inverse UQFT (lUQFT) in e^^ constructed only 
from the positive frequencies of the UQFT in e^^ and had an interpretation from equation (fTSl . The direction of a 
unidirectional signal could be determined directly by the ratio of the Riesz transform components. A scale-localised 
version of this construction would correspond to filtering the image in space and spatial frequency, thus isolating 
structure with a given period and spatial location, and then constructing the Riesz transforms of the local component. 
The resulting quaternionic function could be represented as a quaternionic object. The route to constructing the 
scale localised monogenic signal in one step will be using quaternionic mother wavelet functions. 

Definition 3.3: Monogenic wavelets. 
The monogenic extension of any given real-valued mother wavelet ip^^^x.) is given by 



V'^+Hx) = 7WV'^"^(x) = ^(")(x) +7eV'^"^(x) = ^('■^(x) + i ^)j'^(x) + j V)j"'(x) 



(70) 



where = 7^l V^*^^. = ^2 V'^''^- 
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Note that the terminology 'monogenic wavelet' has been utilized before in a different context: see work by Cerejeiras 
et al. [70]. The monogenic wavelets of [70] are defined on the unit ball and based on Mobius transformations, and 
are distinct from our wavelets. 

Lemma 1: Properties of the wavelet Riesz components. 
Each component of the Riesz transform of a wavelet is also a wavelet. 



Proof: As 



^(e)(f) 



1,2, 



square integrability follows as the wavelets were constructed from a real mother wavelet. Relations (14 U . i42l . 
hold of course, for V'^^^ and its components ip'^^'^ ip^^\ and V'ij^- Denote by c^is) the value of equation d55l) 
with V'/j ^ • Equation (BIT) in combination with the fact that cos^ and sin^ are always between and 1 , imply 
that < c^(i) < c^(r) < oo and < c^(2) < c^(r) < oo. Consequently and satisfy the admissibility 
condition, dSSl . and can be considered as wavelets. The lemma implies that by the monogenic extension of the 
real wavelet we have constructed two additional well-behaved functions and we may then consider the properties 
of the monogenic extension as a whole. 

Proposition 8: Admissibility of the monogenic mother wavelet. 
The monogenic extension of a real wavelet is also a wavelet. 

Proof: It is obvious from dTUll and ^ that |^'(+)(f)|2 = |^'W(f)|2 + |^rg)(f)|2 + |^g)(f)|2 _^ 2 sincp |^'W(f)|2. 
Integrating over all frequencies as in d55t and making use of the Hermitian property of the FT of a real image 
we find that c^(+) = c^(r) + c^(i) + c^(2). Therefore, since each of V'/j^ and satisfy the admissibility 
condition, also satisfies the admissibility condition and can be considered as a monogenic wavelet. Moreover, 
c^(i) + c^(2) = c^(r) , so c^(+) = 2 c^(r) also holds. Hence the monogenic extension of the wavelet is a valid 
analysis wavelet and we may consider its further properties. 

The monogenic wavelets may be represented in polar form via its modulus, orientation, and phase respectively 
given by: 



{+) 



X) = ^[V'^^')(x)]2 + [V.g)^(x)]2 + [^g)^(x)]^ 



(71) 



e^.^(x) 



[V'g,)^(x)]2 + [^g)^(x)]2 

icos(t'g(x)) + jsin(i/5(x)), 



+ j 



[^g)^(x)]2 + [^g)^(x)]^ 



tan 



sgn 



^[V'g,^^(x)]2 + [^g)^(x)]2^ 



(72) 



(73) 



The orientation is thus given an angular representation. The wavelet is a function with spatial energy given by the 
modulus, has a local zero-crossing structure given by (j)^^ (x), and is showing a local orientation preference to angle 
u^^ (x) . The CWT is decomposing the image in terms of localised oscillations with a local period determined from 
0^^(x). For a more thorough discussion of the interpretation of the polar representation of Tp^\-x.), see [71]. 

Proposition 9: FT and QFT of the members monogenic wavelet family. 
The FT and QFT of the translated, dilated, and rotated monogenic wavelet is 



M/(+)(f) 



,(0 



,(1) 



,(2) 



{1+ [-feces I 



+ sm{<))-e)]} 



(74) 



2 



(r), 



[1 + cos(r? -e) + sin(r? - 6)] vI/^-^(r_eq) - fe^'g Vv2+0/2q) 



+ [1 + cos(?? + e)+ sm{7] + 6)] (req) + k^'^> ( J^/a-e/aq) 



'2jiTq2b2 



(75) 
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Proof: See section ITl-AI As noted by Holschneider [32] the entire family of wavelets constructed from a 1-D 
analytic mother wavelet are analytic signals. Analogously, for the monogenic wavelets, we may state the result: 

Theorem 3: Monogenicity of the members of the wavelet family. 
The members of the family of equation dSTl generated from the monogenic wavelet through are also 

^-monogenic signals, this should be contrasted with the 0-hypercomplexing wavelets, that are not 6'-hypercomplex 
signals. 

Proof: This is not a trivial result as we noted in equation (I44t . and corresponds to more than just a simple 
2-D extension of the 1-D result. An arbitrary member of the wavelet family may be recast as V^'^'H^) + 

the Riesz kernels commute with dilations and translations, [VajTZj] = 0, 
[T\^,TZj] = 0. Of course Rg^lj^'''\x) is a ^-monogenic signal/function from equation ( l40b . and thus from proposition 
^ we may note the result. Denote by -0^^^ {x) an isotropic real mother wavelet function. 
I) Isotropic monogenic wavelets: 

Proposition 10: The wavelet family of the isotropic monogenic mother wavelet. 
If = ip^'^^x), then any member of the monogenic wavelet family in the Fourier and Quaternionic Fourier 

domain have a simpUfied form with = [a, 0;b]^by: 



V'^^) (x) = (x) + i [008(0)0-^7^^, + (^)J + [- (^) + cos(0)0^;^ (x 

^W(f) ^ {1+ [-k cos {(f>- 6) + sin {(f>- 9)]} a¥^\af) e-^^""^, (76) 

$W(q) ^ e-^^^^^^^a'f (^)(ag) [1 + cos(0) (cos(??) + sin(r/)) + k sm{9) (- sin(r?) + cos(7?))] e'^''^'^'^^ 

Proof: See section III-BI As an example of the isotropic monogenic wavelet, see Figure |3] The isotropic structure 
of the magnitude is clear, whilst variations over the xi and X2 directions may respectively be determined. 

Proposition 11: 6 = it Equivalence to the anti-monogenic wavelet. 
When the monogenic mother has a radially symmetric real part, then the monogenic wavelet coincides with the 
^ anti-monogenic wavelet, and with = [a, + vr; b]"^ , 

V'^+)(x) =V^-)(x). (77) 

Proof: This we can see by replacing by + vr in (17^ . 

2) Directional monogenic wavelets : The isotropic monogenic wavelets will localise an image to a particular 
scale at each spatial position, and then represent the structure present at that scale and position by a plane wave. 
In some applications it may be suitable to localise in both scale and orientation, as several components are present 
at each position and at the same scale in different directions, and we start by constructing a quaternionic wavelet 
starting with a real directional wavelet defined by equation ( I69l l. The Riesz transforms are found in Appendix III-CI 
and we obtain: 

7^l {0d} (x) = 4 r (5i) (52)1 sin(2^xi(^^)) dh = M^) 

Jo Jo 9 ^ J V2 

7^2 {V-dKx) = 0. 

The complex directional wavelet (see footnote |2l on treating 'tp^j^\x.) as complex-valued) is formed by: 'tp^j^\x.) = 
'i/'D(x) + iV'z),3(x). V'i^'*(x) is a directional monogenic signal. 0;^'*(x)'s directional structure may clearly be 
observed from a contour plot of its two components in Figure 0] ■ip^^\-x.) similar to the complex directional wavelet 
of [1], but is not equivalent to -(/'^^''(x), defined in section llII-A.3l When the Riesz component is constructed from 
'(/'^(x), the spatial frequency modulation to ^'^(f) is fi/f rather than sgn(/i). Of course ^'d (f) is mainly limited 
in frequency to f = [±V2fo, of, and for f = [±V2fo + Sfi, 6/2] , fi/f = ±1 + Sfi/f^ + o{5f^), and 
so up to a small corrective error term V'r>.2(x) = V'd,3(x). This is verified by the visual similarity of Figures ^ 
(b) and (c). The advantage of using a second component constructed from the Riesz transforms, i.e. using 3(x) 
rather than using 2(x), the hypercomplex wavelet second component, is that the quaternionic mother wavelet 
when rotated satisfies the relationships given by equation (l44b . The complex function is then a bona-fide monogenic 
function coupled with a single quaternionic conjugate, an anti-monogenic function. 



',(1) , .^.//3A./.(2) 



',(1) , „..rfl^./,(2), 
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Fig. 2. The real part (a), first, second and third imaginary parts, (b), (c) and (d), of the isotropic hypercomplexing Morse wavelets, with 
n = 0, / = 9 and m — 4. 



C. Quaternionic Morse wavelets 

The 1-D Morse wavelets were introduced by Daubechies and Paul [72], [73] and applied in signal processing 
by Bayram and Baraniuk, as well as Olhede & Walden [33], [74]. In 1-D, according to [33], the analytic 7, n] 
Morse wavelet is '(/'[d „./3 ^(O = '^lUn-P'^i') + J'^S n /3 7(")) and is defined in the frequency domain (c/ [33, 
p. 2666]) by 

*Sn;A,(/) = V2A„;^,,(27r/)^e-(2-/)^L^ [2(2^/)^] , / > 0. (78) 

For / < 0, ^i£)^„./3 ^(/) vanishes. We use the notation c = (2/3 -|- l)/7 — 1, where /3,7,€ K+. Moreover, 
^ra;/3,7 = \/ "^7 2^+"^ T {1 + n) /T {n + c + \) is a normalisation constant and L^(x) is the generalised Laguerre 
polynomial which can be written in a series form: 

Tc. N V , r(l +n + c) 

7^0 ^ r(l + c + r) r(l + n-r) r! ' 

Metikas and Olhede [44] defined isotropic real Morse wavelets, augmented into quaternionic functions, using the 
Riesz transforms. The work in [44] determined specific families of isotropic wavelets that could be considered 
optimally local to a given radial domain, as well as performed estimation of local image characteristics using what 
in this article would be denoted isotropic monogenic wavelets. The n-th order real isotropic Morse wavelet is a 
radially symmetric function, V'i^]rn(')' where / = /3 + 1/2, and m = 7. Its FT and QFT are given by, [44], 

- = (2vr/)' e-(2-/)'" L^'[2 (2vr/)-], (79) 

where c' = (2/ 2)/m - 1, m > 1, / > 0, / > m/2 - 1, and A'^.^^^ = yjm 2'='+i r(l + n)/T{n + d + 1). The 

first two requirements are necessary to make V'i*^lm(/) ^ wavelet, and the second constraint ensures the wavelet is 
related to a suitable localization operator [44]. The 2-D generalization of the analytic Morse wavelets is approached 
by the following method: we first combine the even 1-D Morse wavelets in a suitable tensor product, or we use 
isotropic real wavelet functions as ^(''^(x). We then, to mimic the analytic wavelets, augment the real mother 
wavelet by two or three wavelet functions and form a quaternionic wavelet function, by yVV'^^H^)- 

Definition 3.4: Separable hypercomplexing Morse wavelets. 
The (3, and 7, separable hypercomplexing Morse wavelets are for any (ni,n2) € defined by 

(^l)V'iD,n,;/3,7(^2) " fcV'S,„i ;/3,7(^l )^iA".2 A7(^2) • 

For a plot of the real-valued components of a given separable quaternionic hypercomplexing Morse wavelet see 
Figures ^ (a), (b), (c) and (d). The separable hypercomplexing Morse wavelets exhibit the chequer-board patterns 
of separable wavelets, and can, once the CWT has been calculated, be used to find local separable structure in an 
image at a given orientation. 

Definition 3.5: Isotropic hypercomplexing Morse wavelets. 
The I, and m, isotropic hypercomplexing Morse wavelets are for n G N, defined by 
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Fig. 3. The real part (a), and two imaginary parts, (b) and (c), of tlie isotropic monogenic Morse wavelets, with n = 0, Z = 9 and m = 4. 



For a plot of the hypercomplexing isotropic Morse wavelets see Figure |2 It is clear that these wavelets have an 
isotropic support but recognize variation associated with the xi, X2 and diagonal directions. 

Definition 3.6: Isotropic monogenic Morse wavelets. 
The /, and m, isotropic monogenic Morse wavelets are defined in the frequency domain for n G N by 



(+) 



fccos(,/.)+sin(0))^^':j„(/). 



By using the inverse Hankel transform the real part of the wavelets may be retrieved in the spatial domain: 



/•oo 

/ df f (2^/)' e-^^-/)" Ll[2 (2vr/)"] Jo(2^/x) 
Jo 



The monogenic extension of the n-th order mother Morse wavelet is found by 



V'^l (x) 



{% COS X + J sin x) 27r 



) 

4f / <l^(f) -/iiavr/x) 



(82) 



(83) 



(84) 



Exact forms for these functions for special choices of parameters can be determined analytically (see Appendix III-Dt . 
Figure 15] displays a plot of the isotropic monogenic functions: it is clear that the wavelets have an isotropic support 
but recognize variation associated with the two spatial directions. Note that the Poisson scale space introduced by 
Felsberg and Sommer [47] is not a special case of a isotropic monogenic Morse wavelet as their analysis function 
is not zero at f = 0, and would require setting / = 0. 

Definition 3.7: Directional monogenic Morse wavelets. 
The /?, and 7, directional monogenic Morse wavelets are for n G N defined by 



As expected V^D^o-g 3 corresponds to a complex directional wavelet, and is given by Figures |4] (a) and (b), where 
the other quaternionic components are identically zero. 

Definition 3.8: Hypercomplex Morse wavelets. 
The /?, and 7, hypercomplex Morse wavelets are for any integer valued n defined by 



,n;/3,7 



,'n;/3,'y 



(^2)} 

(X2) - V'S,„;/3,^(a:i)V'S,n;/3,7(^2)}} • 



(85) 



V'S:n;k7(-) 



R 



'7r/4 



R. 



(X2) 

V4 {^S,n;/3,7(^l)'^S,n;/3,7(^2) 
W f^.^ Ao) 



(^2)} 



(86) 



^S,n;/3,7(^l)'^ 



lD,n;/3,7(^2) 



^/j^J^g 3 (x) , is thus calculated from ^1^0-93(^1) ^'^^ ^iD 9 3(^1) equation (l69t . and plotted in terms of its 
real and imaginary component see Figures 0] (a) and (c). The real component is similar to ^^Q.g3(x), but recall 
that the two imaginary components are not equal, i.e. V'D,2,0;9,3(x) 7^ V'r',3,0;9,3(x). Directional wavelets have been 
previously constructed from wavelets related to the 1-D Morse wavelets, see work by Antoine et ah [20], [21], [68], 



but these are distinct from ip^j^l^.f^ ^(x), as well as tp^^^.'p ^(x). The wavelets defined by Antoine et al. [20][p. 319] 
are based on the generalized 1-D Morse wavelets with (3 = I = m and 7 = 1, and are referred to as 2-D Cauchy 
wavelets. In their construction a covariant change of axes has been implemented to make the wavelets compactly 
supported on a cone in the frequency domain, rather than by rotating a function supported in the 45°, or 7r/4 line. 



.,{++) 
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D. Choice of mother wavelet function 

The Morse wavelets may be considered as optimally localised over a given region of space and spatial frequency 
space, see for example Olhede & Walden [33] and Metikas & Olhede [44]. Note that the hyperanalyticizing 2-D 
extensions of the real mother wavelet functions are localised to the same region as the real wavelet functions. 
In order to choose an analysis mother wavelet, we either take a real mother wavelet that is isotropic or one 
which is anisotropic. If the directional choice of localisation is necessary to disentangle local structure, we use a 
combinations of the tensor products of 1-D wavelets. However, if the isotropic choice of localisation is sufficient 
to separate components at given scale and spatial positions, then we base the quaternionic filters on the isotropic 
Morse wavelet. 

For a pre-specified spatial and spatial frequency region of either a directional or isotropic form we choose /3 and 
7, or equivalently / and m, to calculate a suitable family of even or isotropic wavelets that are optimally local to a 
given region (for a more thorough discussion see [33] or [44]). Thus j3, 7, or I and m, should be chosen depending 
on how the image should be localised in space and spatial frequency. As a second stage a suitable quaternionic 
extension of the real-valued mother wavelet is chosen. The choice of quaternionic extension will depend on what 
assumptions can be made regarding the underlying structure of the observed local variation. As may be recalled from 
Figures ^ El IH and |4j each quaternionic extension represents very different local oscillatory structures associated 
with the same period. The separable hypercomplexing wavelets are supported in two orthogonal directions, each 
real component of the quaternionic wavelet capable of extracting even/odd and diagonal structure. The separable 
filters resemble the quaternionic Gabor decomposition suggested by von Biilow and Sommer [49], and the wavelet 
coefficients represent structure separable in a rotated coordinate frame. The directional hyperanalyticizing wavelets 
are constructed from a tensor product of 1-D wavelets, but where the directional wavelets are a sum of 1-D even and 
odd function tensor products, as was suggested by Selesnick et al [1]. Local plane waves, isolated in a given direction 
can be represented well by such a decomposition. The isotropic wavelets mix energy isotropically, i.e. localise an 
image isotropically. Both the hypercomplex and monogenic wavelets are built from the idea that the real wavelet 
isolates a single component by its scale and spatial localisation - subsequently the isolated components will be 
represented by their plane-wave (monogenic), or separable wave (hypercomplex) polar representation. This method 
resembles a Hilbert spectrum representation of a function in 1-D [34]. Once the wavelet transform coefficients 
have been calculated at a single orientation, if the assumption of single isolated structure holds, then no further 
calculations are necessary. This represents a substantial reduction in computational time. Once the CWT coefficients 
have been calculated, the local structure of the analysed image is represented in terms of the quaternionic transform 
coefficients, and their polar representation. 



IV. Quaternionic wavelet coefficients 

Having constructed various quaternionic mother wavelets, each suitable for the analysis of local structure, we 
outline the properties enjoyed by the coefficients. Equation ^ defined the CWT coefficients of image (7(x) analysed 
with a wavelet family constructed from wavelet V'(x). If V'(x) is quaternionic then equation ^ takes the form: 



M^\9) = I I 5(x)V|(x)fi'x (87) 



00 J —00 

00 roo 



/oo roo p 
/ 5(x) V^^Hx) - iijfi^) - - fcvf )(x) 

-00 J —00 



= w^;\t9)-iwf{$;g)-jwli\t9)-kw^^\tg), (88) 

thus defining w^\^;g), u = r, i, j, k. The local structure at ^ is given a four- vector valued representation, via 
some phase/s representations or relative magnitudes of the w^^\$,; g) . 'Local energy' is determined from the local 
magnitude, \w^{£,;g)\ , of the image ^(x). Define C = [«, ^jffe]^ > where f;, is the Fourier variable after the Fourier 
transform in b has been implemented, and k = [a,9,qh\^ , where q^, is the quaternion Fourier variable of b. 

Definition 4.1: Scalogram of g{'x.). 
The local energy of g{x.) at ^ is defined by the scalogram and is given by 

s^itg) = \M^;g)\' = nj^^'itg) + 4^\^;g) + njf\$;g) + wl^'^\tg). (89) 
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Fig. 4. The real part i/'i3.o;9,4(x) ((a)), and first imaginary parts ((b) V'd,3,0;9,4(x) and (c) '!/'-D,2,0;9,4(x)) of the directional monogenic 
Morse wavelet and the hypercomplex Morse wavelet, with n = 0, (3 = 9 and 7 = 4. 



Irrespective of what quaternionic mother wavelet function is used, the magnitude square of the coefficients is 
interpreted as local signal presence. The interpretation is appropriate as the four wavelet functions ip^'^\x.) with 
u = r, i, j, k, are chosen to be local to x = b and scale a, where the localisation of ip^'^^x.) with u = i, j, k, is 
similar to that of ^l)^'^\^x.) from our previous comments in sections HTHlandlirB Note that -V'(°°)(x) enjoys the 
same localization as ^^°°\x.). 

The relative relationships between the other components of the quaternionic wavelet will depend on the choice 
of mother wavelet function used to decompose the image, and so the phase representation of structure will vary 
with the mother wavelets used. To obtain coefficients with interpretable phase relations we calculate g) , so 

that it corresponds to a 0-hyperanalytic signal in index b for any fixed value of 6 and a. The signal g{-) could also 
be represented in terms of some weights attached with functions that are hyperanalytic in x for any fixed 

^0- In the latter case equation ^ states that the signal may be reconstructed in terms of the weighted and rotated 
hyperanalytic functions In 1-D these two perspectives coincide, as analysis of a signal in terms of a analytic 

wavelet yields both analytic wavelet coefficients and a reconstruction in terms of the analytic wavelet functions. 
In 2-D the problem is more complicated, due to the non-zero k component in the hypercomplexing wavelet. We 
introduce the additional notation of Jt^/2I^ = [a, J-K/2^h\^ ■ 



Proposition 12: Forms of the FT and OFT of wavelet coefficients with real wavelet. 

ir] 



The FT and QFT of w^J^\^;g) for a real signal 5f(x) with a real-valued wavelet respectively are 



wi'\Cg) = G(ffe) aM/M*(ar_eff,), (90) 



= ^^<)(-;5) + ^<^(J./2-;5). (91) 

Proof: Equations ( l90b and ( I9U follow by direct calculation, and are only stated so that a comparison can be made 
with the quaternionic wavelet coefficients. 

Proposition 13: Forms of the FT and QFT of wavelet coefficients with quaternionic wavelet. 
The FT and QFT of w^{$,;g) for a real signal g{^) with a quaternion-valued wavelet V'(-) respectively are, 

WAC;g) = G{h) \a^>^'->{ar_eh) - ia^^''^* {ar ^eh) - j a^^^> {ar _eh) - ka^^''^ {ar _gik > 



Ik 1 

W^;Q{i^;a) = -jW^{f^;g) - -W^>.r)_i^(^)j^j^u)-k^w{n]g) + -W^{J„/2i^;g) 

k 

+—W^{r-) _^_j^{f) _i.^(k) ( J^/2k; g) (92) 



/ ^^W^iK-g) + ^W4J^/2K;g). (93) 



— W^^{K-g) + — 

Proof: See appendix IIII-AI Equations ( I92t and (l93l illustrate the new properties afforded by the quaternionic 
decomposition, in comparison to a real-valued decomposition of ^(a;). The i and k coefficients will enable us 
to characterise structure in the xi axis, and we are not constrained to average the coefficients over xi and —xi, 
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cf equation ( 19 U . With a well-chosen quaternionic mother wavelet this will give us tools to represent oscillatory 
structures. 

A. Stability to changes of axes 

The axes of observation are not necessarily aligned with the local axes of variation in the image. The CWT is 
covariant with respect to the transformation x ^ r_6i/x — x', i.e. if we observe g{^) = g2{r^Qi{yL — x')), rather 
than fii2(x), then using equation (2.18) of [21][p. 262] on each of the real valued weights associated with 1, i, j, 
and k, we may formally note that: 

w^{$,; 5(x) = w^{a, e - 0', r^g h - x'; 52(x)). (94) 

Having noted this equivariance of the quaternionic CWT, one might assume that no discussion needs to be provided 
of changes of axes of observation. However, a discretization of the calculation of the CWT, implies that we calculate 
the CWT coefficients at a sampled subset of all values of the locality index ^. Small misalignments in space between 
the wavelet function and the object under observation may cause the CWT coefficient at a given value of ^ to be 
small even if the signal has a large contribution at x = b and f = a^^r^giQ. The down-sampling inherent in most 
filter bank implementations of the CWT exacerbates the spatial initialization problem. In 1-D time shift variance has 
been considered in great detail see for example [22], [24], [75], and small misalignment in time may be considered 
in terms of shifts in phase. In particular, [22] [p. 1794] discusses the relationship between small spatial shifts and 
phase shifts. It is important that, at any given ^, the transform of structure that corresponds to the same space 
and spatial frequency locality should not correspond to very different \w^{^; ^(x))! due to small misalignments in 
space between the image and the wavelet. This becomes equivalent to requiring that the magnitude of the CWT 
does not change under phase shifts of the signal. 

In 1-D it was shown in equation (6.11) and the subsequent discussion in [56] that, if a real signal of sinusoidal 
form is phase shifted by a constant angle Og, the analytic CWT of this signal will correspond to a complex phase 
times the CWT of the original signal. This ensures CWT magnitude invariance to phase shifts, and produces a 
stable CWT. We shall see that the hyperanalytic wavelet coefficients exhibit similar structure, but the definition of 
a phase-shifted signal in 2-D will vary with the choice of local structure. Depending on the local structure of the 
image a phase shift will be defined as a single shift in a given direction, or as two shifts in perpendicular directions. 

B. The hypercomplex CWT 

The hypercomplex signal of a real image (^(x) represents a real-valued image in terms of a four- vector representing 
local structure, where the local variation in the image is considered in terms of the xi and X2 axes separately, as 
well as joint variation in the two axes. When the 0-hypercomplex signal is formed, variations in a rotated frame 
of reference are considered, and this naturally introduces ^_5i(x), rather than g{l^) as the object under analysis. 
If there are separable oscillations present in gf(x) at a spatial position x, in direction along and perpendicular to 
the axes rotated by 6*, then using the polar representation of the 6'-hypercomplex signal given in equation (IIHI 
will allow for the characterisation of the local structure. The local structure is, in this instance, given in terms of 
the amplitude or energy associated with the variation, the period in the first axis' oscillation, given by 5_e^6i(x), 
as well as the period in the second axis' oscillation, given by /?_e e(x). If there truly is a completely separable 
structure present, then 7-6i.e(x) will confirm this assumption by taking the value of zero. If not, a non-zero value 
of 7_e g(x) will characterise the non-separability of the oscillation. By recalculating the decomposition at different 
values of 9 e [0,7r/2], the orientation may be determined from the analysed image, as the value of 6 that yields 
7_e ^(x) = 0. Unfortunately most observed images do not correspond to the simple structure of separable variation 
at any given x, but variation exists at many different spatial scales and orientations at the same spatial position. 
For such images some localisation needs to be implemented when forming the hypercomplex representation of the 
image variation. For an image coiTcsponding to variation at many different scales, the full structure of the image 
at a given spatial position will not be well represented in terms of an amplitude, and two periods. Once different 
structures associated with given variations have been isolated in the image, the structure of each component may 
be given a hypercomplex representation. Under the assumption of local separability the hypercomplex CWT, or 
HCWT, will be calculated. 
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To isolate the individual component we chose to combine the hypercomplex representation of an image with the 
CWT. The wavelet will localise the image in spatial position, scale and potentially orientation, and is assumed to 
be sufficient to separate the different components. In theory, to construct the HCWT, the CWT may be combined 
with the hypercomplex representation in many different orders. We could start by calculating the hypercomplex 
signal annihilating all anti-hypercomplex components in equation dTTt . and then localise the signal in position, 
scale and orientation with a real valued wavelet, with the aim that once localised the CWT would isolate a single 
component that would naturally be given in polar form. We could equivalently argue that we should first project 
the image in position, scale and orientation using a real wavelet to isolate a single component and then find the 
hypercomplex representation of the localised in b signal. It is most convenient to use a single wavelet decomposition 
that produces quaternionic coefficients in one calculation. Fortunately, a single quatemionic wavelet which is not 
itself hypercomplex can be chosen so that both operations can be done in a single step. The three procedures give 
the same filtered image, and this defines a local 0-hypercomplex signal with interpretable magnitude and phase 
structure. To demonstrate the properties of the coefficients it is most easy to work in the Fourier and quaternionic 
Fourier domains. 

Proposition 14: FT of HCWT. 
The CWT of the real signal (7(x) using the hypercomplexing wavelet '(/'^^(•) has a Fourier transform given by 

W^-'^HC, 9) = [l + sgn {[r^efbW] [1 " fegn ([r^effeli)] wj^""\C 9) (95) 

Proof: The result follows by direct calculation. The hypercomplexing wavelet is filtering the signal in the rotated 
coordinate system, constructing a signal locally odd in (r„6ix)i as the i component, a signal locally odd in the 
(r„5ix)2 as the j component and a signal locally odd in (r_5ix)i and (r_0x)2 as the k component. 

Proposition 15: QFT of HCWT . 
The CWT of the real signal (^(x) using the hypercomplexing wavelet ip~^^{-) has QFT given by 

= ^-^(l + sgn([r_eq,]i))(l + sgn([r_,q,]2))PFi^^)(A^;5) 

(1 + sgn ([req;,]i)) (1 + sgn {[rec^b]^)) PF^^'V^s'^; d)- (96) 

Proof: The proof of this calculation follows by direct calculation. Comparing this equation with equation (l25l 
we observe that this is a ^-hypercomplex function, i.e. hypercomplex in a rotated frame of reference in terms of 

F(r^eq) = a^'('^^)*(ar_eq)G'(q), i/(qi) =aM/(^^)(aqi)G(reqi) = avI/(^^)(aqi)G_e(qi). (97) 

The hypercomplex CWT simultaneously constructs a 6'-hypercomplex signal as localising the signal in frequency. 
/i(b) corresponds to the CWT of ^_0(x), calculated at rotation 6* = 0, and we may also interpret w^^^\$,] g) as 
the 0-hypercomplex extension of the real CWT with = of the signal 5f_6i(x), i.e. calculating the real CWT of 
the real signal at kq = [a, 0,b]^ , and then extending it to a 0-hypercomplex signal. By equation (l25l we note: 

«S('^o;5-.))^^^^ = (1 + sgn([r_,q,]J) (1 + sgn ([r.^q^J^)) (i^) wi^'\K-g) 

+ (l + sgn([req,]i))(l + sgn([req,]2)) (^) W^^'Xj^/^n^g). 

Thus w'^^^\$,\g) is a 0-hypercomplex signal and corresponds to constructing the 6'-hypercomplex extension of 

w^^^\$,Q;'g-e{'^))- This gives the interpretation of the coefficients as the appropriate hyperanalytic extension of a 
signal local to = [O)0,/3] existing, and naturally described, in direction 9. 

Proposition 16: The FT of the HCWT of the 0-hypercomplex signal. 
The Fourier transform of the hypercomplex CWT of the hypercomplex signal takes the form: 

VFi++)(C; g'^Xi) = 2 [1 + sgn ([r^efbls)] [1 " ^sgn {[r^eh]^)] wl;'\c, g). (98) 

Proof: See appendix IIII-B I This results is intermediary to proving that the order of operations of constructing the 
hypercomplex signal and calculating the CWT is not important. The appendix gives the form of g^^), 
for Hi and 112 taking any of the values it. For the 1-D signal, only the analytic wavelet transform of the analytic 
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signal was non-null - the analytic wavelet annihilated the anti-analytic component. This is not the case for a 
hypercomplex signal, as in not all three of the anti-hypercomplex components are annihilated. 

Proposition 17: The FT of the real even CWT of the 0-hypercomplex signal. 
The Fourier transform of the hypercomplex CWT of the 6'-hypercomplex signal takes the form: 



(++). 



(1 + sgn {[r^eh]2)) (1 " ^sgn {[r^eh]^)) <''^(C; q) 



(99) 



Proof: See section ITlI-CI We interpret w)p {£,]g_QQ) as the localised 6'-hypercomplex signal of the signal 
observed in a different axes: this should have a polar representation that is representative of the signal structure. 

Theorem 4: Construction of the local 6'-hypercomplex signal. 
By calculating the hypercomplex CWT of a real signal g(x), a local ^-hypercomplex signal is constructed. 

Proof: See appendix IIII-DI and note that 

(++) 1, 



W.. 



(++) 



(++) 



w. 



(++) 



(C;5) 



(++) 



4> 



{++)> 



(ee). 
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(++) 



{++)^ 



<++) 



w. 



(ee) 



(100) 



Once one the above three sets of equations in dlOOt is proved, the rest follow through Fourier, Inverse Fourier, 
Quaternionic Fourier, and Inverse Quaternionic Fourier transforms. Furthermore the equations shed light on the 
interpretation of w^^\^;g). 



'w'^'^\^o'i9-e) is a 0-hypercomplex signal by proposition IT51 thus its polar representation is meaningful, 
as it corresponds to a real signal that has been localised in scale, and extended to a 0-hypercomplex signal. 
Furthermore wll^'^\^;g^_^g^J ) can be viewed as the local contribution of the 0-hypercomplex signal, and is therefore 
also meaningful if there is signal presence at that point. These two quantities are both equal to half w^^^^ g^^g^J ) , 

the hypercomplex CWT of the ^-hypercomplex signal. Thus ■w^^'^\^; g) has an interpretation in terms of 5r_e,6»(x), 
the signal naturally represented in the rotated coordinate system, but observed in another coordinate system. 
w^^~^\^]g) can thus be used to locally at ^ represent a real-valued signal in terms of phase and amplitude. 
Unlike von Biilow and Sommer, we need not assume that a single component is present at x = 5, and it suffices to 
assume that once g{'x.) has been localised to (a, 9) and x = b, then a single component is present. Then w'^^^\$,; g) 
can be represented in polar form using equation (EHl . Furthermore, if we assume that the signal corresponds to a 
locally separable oscillation, the following theorem specifies its representation. 

Theorem 5: The HCWT of a separable oscillatory signal. 
The hypercomplex CWT of an oscillatory signal modelled with x' = r_ex as g{'x) = ag(x) 
cos{2-K(pg^i{x'i))cos{2Tnpg^2{x'2)) whcrc ag(x), Vg,2{-) are assumed to be slowly varying, defining 

^'g,u{'Xi) = l^'fg,u{xi), for ti = 1, 2, and additionally assuming tf'g i{xi) — v?^.2(^2) = C, a constant for all xi 
and X2 such that ag(x) is non-negligible, is given with b' = r_eb by 

w^i^\i-,g) = a,(b)aM/(-) (a/,,i(b,^),a/,,2(b,e)) [cos (2^<^,,i(6;))) cos (27r(^,,2(62)) 
+isin [2-KLpg^i{b[)) cos (271933,2(^2)) + j cos [2'KLpg^i{b'^)) sin (27r(^g,2(&2)) 
+fcsin (27rv9g,i(6;)) sin (27r99g,2(^2))] + o(l) 



w 



(++) 



{i;g) = aa(b)avl/('=^)(a4i(b,0),a42(b,0))e 



27ri</3g,i(fe'i) 27rjV,,2{6J>) 



+ 0(1). 



(101) 



where fg{^,9) = [99^ li^i)-, 9^32(^2)]^ • Note that the constraints imposed on g{^) includes a local plane wave 
structure as then one of the two phase functions will itself be constant, and separable structure with different local 
periods, as this will be modelled via two different local plane waves separated to different values of a and 6. 

Proof: The proof follows by direct calculation, if some care is taken with the conjugation and appropriate 
assumptions are made regarding the variability of the amplitude and phase functions. It is necessary to use a 



STATISTICS SECTION TECHNICAL REPORT: TR-06-02 



31 



hypercomplexing wavelet for the formula to follow. Let fo = arg^^^g max |^'(^'^)(f)| . We note that a simplified 
description of the separable oscillatory signal ^(x) can be determined from the ridge of the CWT cf [43] [p. 394] 
via analysis for values of ^ in a subspace given by |^ : /g(b, 9) = a^^fo} such that the measure of the degree of 
separability of w'^^^^ g) is zero. Direct calculation shows that, when the CWT is calculated at another angle than 
that of the variation in ^(x), the magnitude of the wavelet transform is less than that of equation dlOU . The ridges, 



and thus /^(b, 6), can be determined from maxima in 
fixed a and b varying 



(++)/ 



when the latter magnitude is calculated for 
but to verify the chosen form of the separable orientation is appropriate we additionally 

(b) = 0, where the latter is the 



calculate the polar representation of w']^'^\$,; g) and determine that indeed 7_ 



k angle in the polar representation of w']^~^\$,; g) , see equation (EHl . Once 9', has been determined we calculate 
(b) and f3_Qi ^giih), these two functions then characterising the structure of the locally separable oscillation. 



where 



(++)/ 



w:,: '{i;g) 



corresponds to the magnitude of the oscillation. Signals corresponding to aggregations of several 
separable oscillatory components can be locally analysed as long as the wavelets are sufficiently concentrated in 
space and spatial frequency. For further notes on separating distinct components see [43] [p. 395], where such 
methods mutatis mutandis can be applied in this context. The main focus in this article is not the development of 
additional wavelet ridge methods, but the expression in equation dlOlt shows the structural representation possible 
with the transform coefficients. 



C. Phase shifts in both orthogonal axes 

We define phase shifts in both directions of variation for separable signals. The shifted signal depends on the 
axes of the direction of variation specified via 9, something that will not be the case for the single phase-shift 
defined via the monogenic structure. 

Definition 4.2: Phase shift for a separable structure of orientation 9. 
For a real signal that is separable when viewed in the correct axes of observation, i.e. g{'x.) satisfies equation i32l . 
the signal phase shifted by 6s = (9s,i ^5,2) is defined by: 



1 fys,2) 



cos (27r/?_e,0(x) - 61^,2) 



3f? 



(102) 



where A|^g is denoted the phase shift operator. A|^g shifts the two separable oscillations in cycle by 9s,i and 6*^ 2, 
respectively. 

Proposition 18: Hypercomplex signal of the phase-shifted separable signal. 
For a real signal (^(x) satisfying equation ( I32t the 0-hypercomplex extension of the phase shifted real image in 
direction 9 is given by: 



'~' (++) 
(A0,,e5)_gg (x) 



(103) 



Proof: See section ITlI-EI 

Theorem 6: HCWT of a phase-shifted separable signal. 
The CWT of a phase-shifted signal, where (^(x) satisfies equation i32l . using a hypercomplex wavelet is given by 



w. 



(++), 



Proof: See appendix IIII-FI 

Corollary 1: Magnitude invariance of the HCWT of a separable signal under phase-shifts. 
The magnitude of the CWT of a phase-shifted separable signal is equivalent to the CWT of the original signal. 

Proof: From equation ( I104t we may note that Ag^g?i;(++)(^; g() = e'^^=-^w^^^\^; g)e^^=-^ , and thus it follows that 



(104) 



(++) 



w. 



{++) 



it 9) 



(105) 
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It thus follows that if the original image was separable then the magnitude of the CWT coefficients calculated at 
the correct orientation is stable with respect to phase-shifts. Given the dependence of the analysis on the orientation 
of the signal it is not surprising to find that the analysis strongly depends on the local orientation of the image. Our 
representation of the same image should change substantially if we assume separability in different directions. We 
may deduce that for signals that can be described as locally separable that small changes in the spatial alignment 
between the wavelets and the signal will not cause migration of energy over frequency bands. The result is important, 
as if discrete filter banks are constructed to implement the CWT then ^++(x) phase-shift invariant in magnitude 
implies that the sampled wavelet coefficients will not suffer from substantial spatial shift variance, but there is still 
orientation dependence. 

D. The monogenic CWT 

The interpretability of the hypercomplex decomposition coefficients depends on the assumption of locally separa- 
ble variation present in the image. To be able to retrieve the structural representation of the signal the appropriate axes 
for analysis have to be identified; this in general requires the calculation of the HCWT at the full set of G [0, vr] , 
unless it transpires that the locally preferred orientation is already known. If at a fixed spatial location x = b and 
given set of frequencies there was variation in only one orientation, the structure could be represented by the scale 
localised monogenic signal. If at any given spatial point x, (/(x) corresponded to a single component, then the local 
characteristics of (7(x) could be extracted from g'^^\^)^ and the polar representation of equation ( R^ can be used 
directly to represent the signal. However, in general, g{-x.) will correspond to a multiscale structure at the spatial 
point X, and thus a scale local representation using the monogenic wavelets is necessary to implement in order to 
produce an interpretable polar representation. The monogenic CWT, or MCWT, represents local structure in terms of 
locally unidirectional variation, where the direction of variation can be determined from the transform coefficients at 
a single local point. In general, analysis with a monogenic, but not necessarily isotropic monogenic mother wavelet, 
is implemented i.e. a quaternionic mother wavelet of the form = ip^^\-) = ip^^'^-) + iip^j^\-) + jip^^\-) is 
used. A non-isotropic mother wavelet is a suitable choice as it is not reasonable to assume that the scale localisation 
alone will be sufficient to separate components present in the image. We use an isotropic mother wavelet, when, 
locally, there is only a single unidirectional oscillation present. 

The HCWT were constructed to ensure that the coefficients rather than the wavelets had interpretable polar 
representation. Unfortunately, we found that these two requirements could not be simultaneously achieved. It will 
transpire for the monogenic wavelets that even if the wavelets are themselves ^-monogenic, so are the wavelet 
coefficients, unlike the previous case. This all is due to the zero fc-component in the quaternionic representation. 
The coefficients are represented in terms of an amplitude, an orientation and a local period. Again, to be able to 
establish the properties of the monogenic coefficients, the transform in the Fourier domain is considered in detail. 
We show that the MCWT annihilates the anti-monogenic component, and this ensures that, equivalently, firstly the 
monogenic signal may be constructed and then scale localised; or firstly the image can be scale localised and then 
the monogenic extension of the local signal constructed. Indeed both operations can be implemented in one step 
using the MCWT, and this then establishes the interpretation of the wavelet coefficients. 

Proposition 19: FT & QFT of the MCWT. 
The FT, and the QFT, of the MCWT of a real signal are given by: 

Wl+\C,9) = [l-kcos{(l)b-e) + sm{(l)b-e)]wjl'\c,9) (106) 

^!hQ ('^'^9) = ^^;Q {i^;9) + COS e [cos n + sin Ub] W^^^q (ft; 5) + sin 9 [cos n - sin Ub] («^; 9)j ■ 

Proof: See section IIII-GI From these equation the frequency domain properties of the monogenic wavelet 
coefficients can be determined. Note that from equations (I106t and i42l we can immediately deduce that the 
MCWT coefficients are 6'-monogenic. If the real mother wavelet is radially symmetric, ip^^^-x.) = ^('^^(x), then 
the rotation has no important effect in the definition of CWT, dSTl . This assumption implies that ^^^\f) £ M, and 
a^(^)*(ar_ef) = a'f(^)(a/). 

Corollary 2: FT and QFT of the isotropic MCWT. 



STATISTICS SECTION TECHNICAL REPORT: TR-06-02 



33 



The FT, and the QFT, of the MCWT based on a real isotropic mother wavelet of a real signal is given by: 

WI;\C,9) = ^(^)(aA)G(fb), (107) 

W^^^ (C; g) = [1+ {-k cos(0, -9) + sin(4 - 9))] {ah)G{h) 

= ^(+)(ar_efb)G(f;,) / ^^+>{ar^efb)G{ib) 

= [1+ {-k cos(0b -d) + sin(0b - 6))] Wjf^ (C; g) 

^i^gi'^'^d) = {Ggiq^b) [1 + cos 6* (cos + sinz^f,)] + sin 6* [cos z^;, - sini/fe] iGQ{q)j} a'i?'^''\aqb). 
These relations will enable us to consider the properties of the magnitude of the wavelet coefficients under rotation, 
and we wish to derive certain relations regarding the properties of the CWT in these circumstances. 

Proposition 20: FT of the MCWT 
The FT of the MCWT of a monogenic extension of a real signal takes the form: 

W^f ^(C;5^+)) = [(l±cos(0))GW(ff,)±fcsin(0)G(-)(f5)] av&W* (ar^^fb) • (108) 

The anti-monogenic extension has an equivalent form with the it on the right hand side replaced by =F. 

Proof: See section HiThI and set 9' = 0. 

Proposition 21: CWT of ^-monogenic & anti-monogenic decomposition components. 
The Fourier transform of the MCWT of the extended ^-monogenic and anti-monogenic extension of a real signal 
take the forms: 

w!j;\C;g%\) = {2 + 2sm{cPh-0)-2kcos{cPb-e)} wi'\c,g) 



2GL+Vf,)a*W*(ar_ef; 



= 2WI^\CJ_%), ^^:\U%) = 2w^;\^-7^%) 
WI^\C;±1) = 0, J^\^;±^,) = ^ (109) 
Proof: See section HIl-HI 

Hence the MCWT of a 0-monogenic signal observed in the rotated frame of reference, i.e. g^^Jg{'x.) is twice 

that of the real wavelet transform of the ^-monogenic signal ^_£i e(x). The MCWT of g^_Jg{'x.) is zero. Thus 
the monogenic wavelet annihilates the anti-monogenic component of the image, in perfect analogy with 1-D, 
but we may further generalise the result to consider any combination of monogenic/anti-monogenic wavelet and 
monogenic/anti-monogenic signal. 

Corollary 3: Annihilation using Monogenic & Anti-Monogenic Wavelets. 
The MCWT annihilated the 0-anti-monogenic image and the anti-monogenic CWT annihilated the 0-monogenic 
image: 

wi^\C;g) = lwl^\c-J%), njf\tg) = \wf\t9%). (HO) 

For any other angle than 6 ±1-k,1 Z, w^^^ g) is a linear combination of w^^^ g^^Jg) and w^'^ (^; gg^_Jg). 

Proof: Comparing equation jlll-HI) with equation ( 11061 ) the first result follows directly: furthermore the annihilation 
follows mutatis mutandis from equation (IA-33t . Thus, when analysing a real signal, the MCWT can be used to 
annihilate the anti-monogenic components in the real image, which is in complete analogy to the d = I analytic/anti- 
analytic case, [56], [76]. See equations (6.8) to (6.10), and in particular the unnumbered equation over equation 
(6.10), of [56]. These show that the CWT of an d = 1 real signal with respect to an analytic (anti-analytic) wavelet 
is equal to one half of the CWT of the analytic (anti-analytic) extension of the real signal with respect to the 
analytic (anti-analytic) wavelet. In other words equation (II lOt is a generalisation of the equation over equation 
(6.10) of [56]. This implies that both in d = 1 and d = 2, the CWT of the real signal or image has the same phase 
and one half of the modulus of the CWT of the analytically or monogenically extended signal or image. 

Corollary 4: Annihilation of vr-rotated component. 
If we assume ip^^^ = Tp'^'^\ then with the additional notation of = [a,0 + tt, b]^ , it follows that: 

<^^(C.;5) = \ w^^\c.J3,o)^ ^f\^.;9) = \ n,f\i.;g%)- (in) 
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If an isotropic wavelet is used then the monogenic wavelet rotated by vr is identical to the anti-monogenic wavelet, 
see equation illl . and so the result is an immediate consequence of the previous proposition. 

Theorem 7: Construction of the local ^-monogenic signal. 
The MCWT of a real image g(x) is equivalent to the scale localisation of the ^-monogenic extension of the signal 
and the ^-monogenic extension of the scale-localised version of the image (/(x) : 

(+) 



Proof: See appendix IIII-II Thus all three local descriptions, i.e. the MCWT of a real signal, the real CWT of a 
6'-monogenic signal and the monogenic extension of a real CWT of a real signal of a real-valued signal may all be 
viewed as equivalent. The left hand side of equation jl \2\ simply corresponds to the MCWT of the real signal, whilst 
the middle equation is the ^-monogenic extension of the scale local signal ^-^(x), which corresponds to adding 
some suitable components to the real part of the CWT of ^(x). The right hand side of equation (II 12t corresponds 
to first forming the monogenic extension of the signal, and then scale-localising the extension to make the phase 
description of the signal interpre table. Hence the MCWT can, assuming a single component has been retrieved at 
local index ^, be represented in polar form using equation ( I46l l where the magnitude, phase and orientation are 
interpretable in terms of a local univariate variations and the ^-monogenic representation. Furthermore if we may 
assume the signal corresponds to a local plane wave the following theorem specifies its representation. 

Theorem 8: The MCWT of a plane oscillatory signal. 
The MCWT of a single component separable oscillatory signal modelled by 

^(x) = ag(x) cos (27r99g(x)) where ag(x), as well as (/'g(x) are assumed to be slowly varying, is given by: 

= ag(b)a ^W(ar_0/g(b)) e2^^"<'('''«>[^«(^)-^^'('^^-«-^«(''))] +o(l), (113) 



where /^(x) = V99g(x) = /^(x) [cos((/)g(x)) sin(</>g(x))] 



T 



n 



e) = [cos(0g(x) - 6) sin(0g(x) - e)Y , 1'('')(f) = l^'^^'^f)! e-^'^^^'^^^), and note that Hermitian symmetry 
imposes ll-^^'^f)! = I'^'^'^^-f)! as well as v9^(f) = -(^^(-f). 

Proof: This result was shown for the special case of the isotropic Morse wavelets in Metikas & Olhede [44], 
and follows by direct calculation. Let fo = argj^^o l^'^'^^f)! • We note that a simplified description of 
the oscillatory signal g{'x.) can be determined from the ridge of the CWT cf [43] [p. 394] via the subspace of 
the locality index given by : fg{h) = a^^rgfo}. Equation (I113t can be used to characterise the oscillation at 
X = b and if the signal would be more appropriately modelled as an aggregation of of oscillatory signals as long 
as they are sufficiently separated, i.e. as long as the wavelet is sufficiently narrow in space and spatial frequency to 
separate the different component. Further analysis of the signal such as that of Olhede and Metikas [44], [71], [77], 
using isotropic monogenic wavelets, can then be implemented for arbitrary monogenic wavelets. Using directional 
monogenic wavelets will allow us to analyse a larger class of signal as it is more reasonable to assume an image 
corresponds to a collection of plane-waves if we localise in scale, position and orientation. 



E. Phase shifts for Unidirectional Variation 

Section IIV-AI noted that the stability of the CWT coefficients under given affine transformations of the argument 
was important. For the HCWT we remarked that the coefficients exhibited desirable properties for phase-shifted 
separable structure. With the monogenic transform we can obtain desirable properties for an arbitrary phase-shifted 
signal, and in some instances obtain invariance under rotation. As the local structure is described as a plane wave 
with orientation determined from the signal under observation this enables us to derive the additional results. 
Furthermore the phase-shift is not dependent on the choice of axes: the monogenic wavelet is identifying the 
direction of variation from the observed image, and then the phase-shift is given in terms of this direction. This is 
not equivalent to the hypercomplex wavelet coefficient structure. 

Proposition 22: Local magnitude invariance under rotation. 
If the real mother wavelet is isotropic, then the MCWT coefficients of a real image g{-K) at ^ will have a magnitude 
invariant to the value of 6. 
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Proof: See section UlLJl Hence the magnitude of the wavelet coefficients does not depend on the rotation angle, 
and we may deduce the result. This implies that we are not sensitive to loss of magnitude due to local orientation 
misalignment. However, we would naturally like to also be able to retrieve the local orientation, from the MCWT. 

Definition 4.3: A directional signal. 
A real signal ^(x) is considered as one-dimensionally directional at frequency f if 3(j)o G [0, vr) such that 



G(f) 



Gif) 



(114) 

G(-f), and thus 



We shall in this article be analysing real images: as g{x) e M its FT is Hermitian G*(f) 
necessitates having both a delta distribution component at (pQ and (pQ + tt. 

Theorem 9: Directional selectivity of the MCWT. 
If the real mother wavelet is isotropic, and the analysed signal g{-x) is directional over the support of the monogenic 
wavelet at ^ with directionality ^o, then with = [^^ 4>o,t>f" 



w, 
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Proof: See section IIII-KI Thus with an isotropic mother wavelet, the magnitude of the Riesz component wavelet 
coefficients is invariant to rotations, but is maximised in the first component if we rotate the wavelet to align 
with the directionality of the variations. Hence we may determine the local directionality of the signal from the 
monogenic isotropic wavelet coefficients. 

1) Phase-shifted plane wave signals : For a plane wave present at a given point x we define the phase-shift 
operation as follows. 

Definition 4.4: Phase-shift of a plane wave. 
For a real signal g{x) = ^_g(r_gx) we note from equation ( l45t that the image may for any rotation angle 9 be 
written in sinusoidal form and the phase shifted by Og version of the signal is defined via: 



W(x) cos (27r0(x) - 0,) 



,{x)e.„(+) 
9-e.i 



X +e= 



3f? {b^+Hx)! exp ((27r<A(x) - 0,)ep_,_, (x)) } 



9 



(-) 



X 
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where Ag, is denoted the phase shift operator. 

Proposition 23: The ^-monogenic phase-shifted signal. 
For a real signal (^(x) the 6'-monogenic extension of the phase shifted real image if zJ^g g(x) is varying sufficiently 



slowly, i.e. assuming v^0^g{-x.) 
form: 

(+)/ 



j^_e)^6»(xo) for all x such that |(7^_^^^(x)| is of non-negligible magnitude, takes the 



iH9)ri^) 



i^iy;(x) 



(x) 



COS (27r0_6»,6i(x) - 9s) + e^_o,o (x) sin (27r0_6 

giV.ei^). 



(117) 



If er. „ „ (x) — ej)^^, i.e. the orientation is constant across the width of the image, then as the Fourier transform 
corresponds to a right-hand multiplication, we may note 

T[{AeJ_%ix)} = e-^^^^— ^^^GL+^f). (118) 

Theorem 10: MCWT of a phase-shifted signal. 
The MCWT of a phase-shifted signal, where the orientation of the signal is constant over the width of the wavelet 
is given by 

wf\t^e.9) = e'^^^-^^o^''^'^wf\t9)- (119) 

Proof: See section IIII-LI This result implies that the phase shift between two images that locally correspond to 
plane waves may easily be determined by equation ( II 19t . 

Corollary 5: Magnitude invariance of MCWT under phase-shifts. 
The magnitude of the MCWT of a phase-shifted signal, when the local orientation of the signal is stable over the 
width of the wavelet is equivalent to that of the not phase-shifted signal. 
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Proof: From equation dl 191 ) we may note that wi^j (^Sff) = e'^'"^"e^^)^^w^^^\$,; g), and thus it follows that 



(±) 



(120) 



Hence the magnitude of the MCWT is invariant to shifts of phase. This implies that there will be no migration 
of energy across scales due to misalignment between the wavelet and the signal. This shows the stability of the 
MCWT to misalignment between the observed signal and the analysis filters. 



V. Conclusions 

In this paper we have defined and provided suitable 2-D extensions of 1-D analytic wavelet decompositions that 
we denoted hyperanalytic. We have constructed classes of 2-D mother wavelets, as well as specific examples of 
such functions that may be used to calculate hyperanalytic decomposition coefficients, and discussed properties 
of the hyperanalytic wavelets and wavelet coefficients in detail. We have stressed the importance of magnitude 
invariance of the coefficients to phase-shifts of the image, and have constructed wavelets yielding coefficients that 
may be thought of as 'locally hyperanalytic' to exhibit such invariance. 

We argued that the design of 2-D hyperanalytic decompositions in continuous space must start by understanding 
2-D hyperanalytic signals. Hyperanalytic signals are constructed as the appropriate limit of hyperanalytic functions. 
A hyperanalytic function is a set of functions satisfying any generalisation of the Cauchy-Riemann system. When 
considering analysis in dimensions higher than 1-D there is in general more than one such system, and so a choice 
of the hyperanalytic system must be made. We considered two different systems in 2-D, namely the monogenic and 
hypercomplex systems. We carefully distinguished between hyperanalytic signals (limits of hyperanalytic functions) 
and hyperanalytic functions in our discussion. 

Once a choice of hyperanalytic has been made, the hyperanalytic signal had to be suitably merged into the CWT 
to yield a hyperanalytic CWT set of coefficients. Two additional problems of interest were then resolved. Firstly the 
rotation operator does not commute with the convolutions that combine to construct a hyperanalytic signal from a 
real-valued image. To determine the structure of rotated hyperanalytic objects in a suitable frame of reference, we 
defined the notion of ^-hyperanalytic as a function that is a hyperanalytic signal in a rotated frame. This definition 
allows for the interpretation of the hyperanalytic wavelet coefficients in terms of local variational structure in 
a given orientation. Furthermore the monogenic and hypercomplex wavelet transforms were related to the local 
UQFT (a unidirectional oscillatory representation), as well as the local QFT (a separable oscillatory representation), 
of an image. We demonstrated that our choices of hyperanalytic wavelet decompositions thus correspond to the 
representation of local image structure as either a plane wave (locally 1-D structure) and a local version of the 
UQFT basis elements, or as a combination of locally separable structures with a natural axis of observation and a 
local version of the QFT basis elements. 

Given a hyperanalytic signal is a three or four vector-valued object it admits a representation as a quaternionic 
function. The local energy of the real signal was described by the magnitude of the quaternion, whilst its local 
structure was described by some suitable phase function/s defined in terms of the quaternionic elements, thus 
parameterising the local structure of the image in terms of local magnitude, scales and orientation of variation. 
Different polar descriptions were chosen for each of the two 6'-hyperanalytic constructions, as each hyperanalytic 
definition is suited to particular forms of local signal variations, and the interpretation of the magnitude and phases 
were discussed. 

Secondly, as we demonstrated, a decomposition constructed using a hyperanalytic wavelet does not in general, 
even for zero rotation of the mother wavelet, yield wavelet coefficients that are hyperanalytic extensions of real- 
valued wavelet coefficients. In fact the wavelet coefficients need not be the real wavelet decomposition of a signal 
that has been extended to a hyperanalytic object. We introduced the concept of an hyperanalyticizing wavelet, 
as a function constructing hyperanalytic wavelet coefficients. In 1-D the analytic wavelets localise the signal in 
scale and time simultaneously as yielding a local interpretable description of the signal in terms of the local 
amplitude and the phase of the signal at that scale. As both the operation of localisation and construction of an 
analytic object commute in 1-D (both are multiplications in terms of complex objects) this simultaneous operation 
is not problematic. This, we demonstrated, is not the case in 2-D (quaternionic objects in general do not commute 
and, in addition, the rotation and the construction of hyperanalytic images do not commute). We showed that the 
hyperanalytic coefficients were equivalent to either constructing the real wavelet decomposition of the real image and 
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then calculating the ^-hyperanalytic extension of the wavelet coefficients or to constructing the 0-hyperanalytic signal 
and then constructing the real wavelet coefficients of this object. This gave the interpretation of the hyperanalytic 
wavelet coefficients, as the ^-hyperanalytic signal of the local projection or equivalently the local projection of 
the 0-hyperanalytic component. To be able to prove these properties we defined the decomposition of any real- 
valued image into ^?-hyperanalytic and ^-anti-hyperanalytic objects, in direct analogue with the analytic/anti-analytic 
decomposition in 1-D [41]. 

We demonstrated that the hyperanalyticizing wavelet either eliminates, or converts, the anti-hyperanalytic objects 
into hyperanalytic objects. The hyperanalytic decomposition allows for the determination of the CWT coefficients of 
phase-shifted images, where the elimination/conversion of the anti-components is pivotal for easy parameterisation of 
the coefficients of the phase shifted signal in terms of the coefficients of the original signal. For both the monogenic 
and hypercomplex coefficients we retrieved magnitude invariance of the wavelet coefficients under phase-shifts, even 
if in the hypercomplex case some restrictions had to be placed on the observed signal. Magnitude invariance to 
phase shift is an important property when discretizing the implementation: the local representation will be sampled 
and so small shifts in space should not change the format of the representation substantially. In the monogenic 
case we additionally demonstrated local invariance of magnitude under rotations with certain classes of mother 
wavelets and proved directional selectivity of the quatemionic coefficients. The stabiUty and other properties of the 
transform depend on constructing wavelet coefficients that are hyperanalytic: this justifies our careful development 
of the hyperanalytizing wavelets. 

In 1-D, analytic wavelets are also used for the decomposition and characterisation of non-stationary oscillatory 
signals, using wavelet ridge techniques [4], [5]. Usage of a localisation method is necessary for the analysis of 
multi-component non-stationary oscillatory signals, as for such signals the phase and amplitude description of the 
local structure in the spatial domain is not informative [30]. In 2-D we showed that wavelet ridges for unidirectional 
and separable oscillatory signals can be determined from the monogenic and hypercomplex wavelet coefficients. The 
advantage of using quatemionic wavelets is that separable oscillations can easily be parameterised in both directions 
simultaneously, and the analysis of unidirectional oscillations is simplified as once components have been separated 
in phase space the orientation of the variation can be determined from the monogenic wavelet coefficients without 
calculating the transform at all values of 6. Our work thus complements existing ridge analysis using complex 
wavelets in 2-D. Given all the components of the quatemionic wavelet can be argued to have the same spatial and 
spatial frequency localisation only a few of the wavelet coefficients of most deterministic image features will be 
of non-negligible magnitude, whilst most noise spreads across all coefficients. Thus the localisation also increases 
the signal-to-noise ratio of the signal in the coefficients that contain most of the signal presence, and this further 
facilitates analysis of noisy realisations. 

Finally, even if general results are of interest, to be able to implement the CWT, explicit examples of mother 
wavelets must be provided for analysis. Existing discrete hyperanalyticizing wavelets [2] and [1] were given an 
interpretation in terms of the general framework developed in this paper. Furthermore we provided families of 
continuous wavelets that are hyperanalyticizing. We chose to base the wavelets on the 1-D and 2-D generalized 
Morse wavelets, that can be shown to enjoy optimal localisation in phase space: we provided five classes of wavelets, 
the hypercomplex Morse wavelet, the separable and isotropic hypercomplexing Morse wavelets, as well as the the 
directional and isotropic monogenic Morse wavelets. 

The study and development of analytic 1-D filters have lead to the synthesis of a wealth of signal processing 
methodology and applications. Before the discrete implementation was developed, continuous analytic wavelet 
decompositions using analytic wavelets were introduced and championed, such as wavelet ridge methods, and local 
decompositions of multivariate time series. This paper demonstrates the potential for hyperanalytic (quatemionic) 
methods, and the representation of local 2-D structure in terms of quatemionic objects. It is anticipated that methods 
based on quatemionic wavelet coefficients hold many novel challenges and potential future developments. 
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Appendix I 
Hyperanalytic Properties 

A. QFT of g^^^\i^) 

Note that G(f) commutes with j as it is a function of j only. We obtain: 

= ^-^G(r_,q) + i±^G(r_, J./^q) + i (^cg^ (r_,q) 

+4^G2)(r_,J./,q)) + (i^Gg)(._,q) + i±^Gg)(r_,J./,q: 

+i (i^Gg)(r_,q) + i±^Gg)(r.,J./,q)) 
1 — fe 

= (1 + sgn ([r_0q]i)) (1 + sgn i\r^ei^^) G(r_eq) 

+ (1 + sgn ([req]i)) (1 + sgn ([req]2)) G(r_9 J,/2q) (A-1) 

Furthermore, if we had started by defining x' = r_ex and then found the QFT in x', with QFT variable q' we 
would have obtained Tq j^e^^''} (^') = (1 + ^gn ([q'Ii)) (1 + sgn ([q']2)) G'Q(q'), and so in the rotated frame of 
reference the QFT is supported on positive quaternionic frequency only. 

B. The 6-Monogenic Signal 

With p(x, y) and qs{^, y) defined by equation ( I34t let us denote by operator notation the construction of a system 
defined from ^(x) satisfying the Riesz system of equations: Vg{'x.,y) = y) * (x), and Qsg{'x.,y) = 

{qs{-,y) * (x). Whennij„g(x,y) = Vg{^,y), and i;S^^i^(x, y) = Qsg{^,y), we denote A; J.^ (x, y) = UR,g(x,y)+ 

i v^^^g{x,y) + j u^]^(x,y) = V ^(x) + i Qi c/(x) + j Q2 g{^), which for y ^ 0+ in the upper half- 
space reduces to the monogenic function y(+)(x) = y(x) + i <^^^(x) + j y^^(x). Rotate the solution to form, 
Rg V g{^,y) + i Rg Qi g{^,y) + j Rg Q2 g{^,y)- Rotations require some attention; they commute with the 
convolution with the Poisson kernel [fig ,V] =0, but not with operation of the convolution with the conjugate 
Poisson kernels [Rg, Qj] ^ 0, see [38]. In particular, 

Re Qi 5(x,y) \ (0,1 Rgg{y.,y) \ / Qi g{r^gx,y) 

Rg Q2 5(x,y) J ' V 22 Re g{^,y) J ' \ Q2 g{r-gx,y) 

As noted by [78], the Riesz system is invariant under dilation, translation, and simultaneous rotation of (xi,X2) 
and {Qig, Q2 g) by the same angle. In other words, (IA-2t is also a solution of the Riesz system, but this time 
defined in variable x' = r„5ix. For y ^ 0+, (IA-21) reduces to 

Rg g{'x) + i [cos 6 IZi Rg g{x.) + sin^ 7^2 Re 5(x)] 

+jf [-sin 6* 7^l Rg g{x) + cos 6* 7^2 Re 5(x)] . (A-3) 

This can be recast as Rg g{-x.) + i Rg TZi g{'x.) + j Rg 7^2 5(x), see [38] for the algebraic properties of the Riesz 
kernels, and in analogy to iA-2l . 

Ug 7^l g{-K) \ ^ ( UiUg g{-K) \ ^ ( IZi g{r^gx) 
Kg 7^2 <7(x) ) \n2 Ug g{x) ) \n2 g{r-gx) 



(A-2) 



(A-4) 



STATISTICS SECTION TECHNICAL REPORT: TR-06-02 



41 



Therefore, if a given function ^(^^(x) is a monogenic signal, i.e. the Umit of a monogenic function, then the 
function generated by rotating its axis by 6 is also a monogenic signal, i.e. it is the limit as y — > 0^ of a system 
of functions, satisfying the Riesz system in a rotated frame of reference. A key component in the functionality of 
the set-up is that |r_ex| = |x| . Hence the 6'-monogenic signal is aptly named. Mutatis mutandis the proofs extend 
to the anti-monogenic case. 



Appendix II 
Properties of Quaternionic Wavelets 

A. FT and QFT of the Monogenic Wavelets 

We note that 



(+) 



ae 



-2i-Kiqibi 



.FQi?,{v^M+#g)+Vg^j}(aq) 



as each real wavelet is transformed according to: 

/oo poo -| 
'OO J —oo ^ 



oo J — oo 
oo roo 



r-e 



x-b 



g -2TTi{qibi+aqiyi) ^^{r) (^_gy) Q-'2-^3{Q2b2+q2ay2) ^2^ 



oo J —oo 



Note that from [79] for V(x) G 



^Q(q) = ^ir^^(q) + ^^^((-^1,92)), 



and we shall use this property on numerous occasions. We find with ^^^^(x) = Rq {ip^^''} (x) that 

.FQi2,{v^M}(q)=.FQ{V^W}(q) 2 
We also note that ^'^'^■'(q) = ^'^''^ (r_0q). Thus we find 



-%'{ci) + —-¥g>{J^/2Cl}. 



2 

1 - k 



+- 



2 

1 + k 



2 '^Qir-9Jn/2<l) + ^-^'^Q\jn/2r-9J^/2CL) 



l + fc^(r) 

k 



^4^V-.q) - ^^gV./2+./2q) + ^^gVv2-./2q) 



(A-5) 



(A-6) 



(A-7) 



Similarly to ^'^'^g(q) we can now derive the QFTs of the members of the families produced from the first and 

second Riesz components of the monogenic mother wavelet, ^/j^^Q(q) and ^^^^^(q). as these two functions are 
real, see also [79]. As before, the crucial step is to calculate the QFT of the rotated real wavelets, by combining 
equation iA-7l with equation (l37l . This yields: 



cos(7? - 9) ^g^(r_eq) - k^'^> {J^/2+e/2^) 



{r), 



+ (-^)cos(7? + 0) fc^g)(J^/2-e/2q) + ^^'^M 



(A-8) 
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1 

2 
1 

+ 2 



^gV-eq) - k^'Q'{J./2+9/2Cl}\ (-J)sm(r? - 6) 



(A-9) 



Putting equations (IA-6t . iA-7t . (IA-8t . and ( IA-9t together, we get the QFT of a member of the family generated 
through translations, dilations, and rotations from the entire mother monogenic wavelet: 



^g-27rigifei 

2 



(1 + cos(?? -e) + sin(7? - e)) f ^g'V-eq) - fc^^'(J./2+e/2q: 



+ (1 + cos(r/ + 0) + sin(?? + 9)) ( ^'[^"^(req) + fc^^'^(J^/2-e/2q 



-2TTjq2b2 



(A-10) 



B. r/je Wavelet Family of the Isotropic Monogenic Mother Wavelet 

From equation i74l we may note that the FT of a translated, dilated, and rotated monogenic wavelet which 
originates in a monogenic mother wavelet with radially symmetric real part respectively is given by: 



{1+ [-fccosi 



+ sin (0-0)]} a¥'\af) g-^'^-f'^ 



(A-11) 



/oo roc 
/ {1+ i-i j cos {(l)-e)+ sin {(j)- 6)]} a¥'\af) e-^^''^"^^-''U^jc 



-OO J — oo 

(x) + i [cos(0)V'5J,^^^(x) + sin(0)V'g)^^(x) 



+ j 



sin(0)V^«^(x)+cos(0)V^)j^;^^(x) 



(2) 



(A- 12) 



This gives the stated results. For the QFT expression, using equation ( I75l l. we note from equation (lA-lOt . recalling 
that the QFT of an isotropic real wavelet is a real-valued object: 



(+) 
Q,4 



[1 + cos(7? -6)+ sm{7] - 6)) ( ^^Q\q) - k^'o'{q) 



+ (1 + cos(r? + 6)+ sin(?? + 9)) + k^'q' {q) ^ 

= ae-'^^^'^'^'^^o'iq) (1 + cos(e) (cos(r/) + sin(7?)) + ksm{e) (- sin(r?) + cos(r/))) e-'^^^'^^^^ 
This yields the desired result. 



-27rj<j2fe2 



C. Directional Monogenic Wavelets 
We recall that 

^^Uf) = -j^^Uf). if />o, ^.s°2(/) = j^s^2(/)> if />o. 

Direct calculation yields the result 

/: 



(A-13) 



{^d} (x) 



oo fOO 



' — oo J —oo 3 f 
^g2,Vf-x ^2f 



V2 



V2 



V2 



V2 



oo /•oo 



^0 

oo roo 



Jo 



91 +92 

V2gi 
9 



^iD (51 ) *S (52)J sin(27rg'^r,/4x) 
^^il i9i) (52)1 sin(2vrxi(^l±^)) d^g. 



V2 



(A-14) 



Thus TZi {iPd} (x) is constant in X2 and odd in xi, and a highly directional wavelet. We can in a similar fashion 
calculate the second Riesz transform to retrieve. 

1 



7^2 {^d} (x) 



V2 



(V4,2) , 



0. 
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Furthermore note that: 



V'd(x) 



OO POO 



OO J — OO 



(e) f fl-f2\ ^(e) ( h+f2 



ID 



V2 



ID 



V2 



(o) f fl- f2 \ ^{o) [ h + f2 



ID 



V2 



ID 



V2 



JO 



^Sd (5i) ^Sd (52) cos(27rg^r,/4x) 



D. Exact form of the spatial representation of the isotropic Morse Wavelets 

To be able to calculate the spatial domain representation of ijj^^'^ (x) , given by equation (l84l l in the instance of 
the Morse wavelets, we make use of the following integrals, see [80], 



2tt 



smm(j) = 2TTj"^Jm{z) s'mmO, 



2tt 



^jzcos{9-<p) ^.os^0 ^ 2^j'"J„(z) COS m0. 



In special cases n = 0, / > —2, m = 1,2 the analytic forms corresponding to equations ( I83l l and ( 1841 can be found. 
Note that Lq(x) = 1. For m = I and n = we get 



1 



r(/ + 2) 



1 



2vr [l + x2] 

where we have used equation 6.621.1 of [81] and Pu are the Legendre functions of degree v and order p, see e.g 
[82]. The monogenic wavelet is 



^o^iW = V'o;/,i(a^)+ cosx + j sinx) 



27r 



For m = 2 we get: 



[1 + x^] 2 Vl + 



where we have used equation 6.631 of [81]. For I = 2, s = 1/2, this wavelet is identical with the mexican hat, see 
e.g. [21]. The monogenic wavelet is: 



A' 

4^,2 W = (* cosx + j sinx)^^ 



/ + 3 



X -r(— ) iFi( 
4 ^ 2 ^ ^ 2 



;2; 



Appendix III 

Properties of the Quaternionic Wavelet Coefficients 

A. QFT of the wavelet coefficients with an arbitrary quaternionic wavelet 

The CWT with respect to a quaternionic wavelet is written in terms of its real components and calculate the 
QFT component by component as follows: 

•^Q W(^;<?)} = {w'^\i;g)}-iTQ {wl\i;g)}-TQ [wf {i; g)} j - i^Q J 

Ik 1 

= -W^{K;g) - -W^(.)_i^(i)j^j^u)-ki,(''){i^]g) + -jW^{JT,i2n;g) 



1 — k 1 + fc 

/ -^W^{i^]g) + ——W^{J^I2n]g). 



(A-15) 
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B. The Fourier Transform of the Hypercomplex wavelet transform of the 6-Hypercomplex signal 

To establish the properties of the wavelet transform, we find the Fourier decomposition of the real signal in terms 
of the 0-hypercomplex and anti-hypercomplex components. We note that for /ii = it, and /U2 = ±, 



w, 



Of course for any real-valued function /i(-) by the linearity of the wavelet transform: 

^Ip^'^'h^; h) = w^;'\t h) - f,iiw^;'\t h) - fi2jwl;''\t h) + ^x^^X2kw'~;''\i; h). 

Thus 



VF;++^(C; h) = WI;'\C h) - iW^;'\C; h)jsgn ([r_,f,]i) - jWf\C, h)jsgn {[r^ehW 
+kW^''\C, h){-l)sgn {[r^eibW sgn {[r^ehW 

„ „-P. 1 Ac^i-^/'U „f, 1 M TI7 



= [1 + sgn([r_ef6]2) " ksgn {[r ^eh\i) - fcsgn ([r„efb]i) sgn ([r^effels)] W^^"\C;h) 
= [l + sgn([r_ef6]2)][l-fesgn([r_ef6]i)]W^i''^(C;/i)- (A-16) 
Thus the FT of the CWT of the hypercomplex image with respect to the hypercomplexing wavelet is given by: 



If we consider signals that are naturally aligned with the wavelet ^/J'^^{■), i.e. take 6' = 6, then a tedious but 
straightforward calculation yields 

wl-'+\C,g'LtJ) = 2[l + sgn{[r,eib]2)][^-ksgn{[r^efb]i)]W;;'\Cg). (A-18) 

This quantity is non-zero in the first quadrant of the QFT domain, when 6 = 0. The hypercomplex wavelet 
transforms of the other three single quadrant supported signals follow from similar computations and are given by 

W^'-+\cJSe!J) = 2[l + sgn([r_effc]2)][l + fesgn([r_,ffe],)]M/^^^)(C,<7), (A-19) 
wi;^\C;g%^) = 2[l + sgn([r_effe]2)][l-fesgn([r_ef6]i)]T^^''^(C,5), (A-20) 

^^^"'''^(C;^^) = -2[l + sgn{[r_eib]2)][l + ksgn{[r.eitW]wl''\c,g). (A-21) 

Thus we have considered the hypercomplex wavelet transform of the ^-monogenic decomposition of the real-valued 
signal fi((x), 

C. The Fourier Transform of the Even Wavelet Transform of the 9 -Hypercomplex Signal, and Coefficients Equiva- 
lences 

One possible strategy to construct an interpretable split into amplitude and phase descriptions would be to first 
construct the 0-hypercomplex extension of ^_0(x) and then take the real wavelet of this object to project the 6 
hypercomplex extension into a small region of space and spatial frequency. The object of this procedure would be 
to extend the aggregate of sinusoids into a suitable quaternionic extension, and then once the quaternionic extension 
has been projected, the projected components are well represented in polar form. It follows that: 

4'^ it g^Xe ) = 9-0,) + iwi;'^ it g%j,) + g%.,H) + 9%.,h) 

+kwl''\t 5)(-l)sgn ([r_effe]i) sgn {[r.efbW 
= {l+sgn{[r_eib]2)){^-ksgn{[r_efb]i))wl''\c,g). (A-22) 
We used the expression of the Fourier transform of a 0-hypercomplex signal given by equation i24l . 
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D. Construction of the Local 9-Hypercomplex Signal 



Firstly let us derive an expression for 
that: 



(++) 



to show equivalence to the other expressions. Note 



oo roo 



(A-23) 



oo J — oo 



Thus if we take the 6'-hypercomplex extension of w^'^\^o, g-e) we obtain the representation: 



(ee)/>. ~ \ 

'{$o,g-e) 



(++) 



wii, '{io,g- 



+3 



w\ '{^o,g^e) 



+ 1 

(2) 



w\i, '{^o,g-e) 



+ k 



(1) 



iee) f J. ~ N 



(3) 



(A-24) 



We thus using equation (I24t we retrieve the following in the Fourier domain: 



iee) f ~ ^ 

w\ '{io,g-e, 



[1 + sgn((r_,f,)2)] [1 - ksgn{{r_efb)i)] G%)a¥^'\ar_efi 



and thus we may note from equations ( l95t . (IA-22t and finally (IA-25t that 



w 



{++) 



(ee)/^ ~ N 



{++) 



(A-25) 



(A-26) 



Finally as the wavelet transform is linear we may note that for any fixed value of the rotation parameter 9 = 9', 
that by equation d27t : g{-x.) = J2±lf'-e^l'i'^)^ we retrieve: 



w 



-(-+) 



{++), 



-{+-) 



,(++)/ 



) 



Hence we find that using equations (IA-18t . (IA-19I I. (IA-20t and (IA-2U that the hypercomplex wavelet transform of 
5((x) viewed at 9' = 6, via equation (ETl the wavelet transform of g{x.) is decomposed into the transform of each 
component, but that naturally sum to a hypercomplex signal: 

W^-^^\C,g) = [l + sgn{[r_eib],)][l-ksgn{[r^g{bW]W;^''\C,ge). (A-27) 

Thus equation (IA-271) naturally agrees with equation ( l95t but by carrying out the calculation we may note that 
we annihilate the sum of the anti-hypercomplex components corresponding to xi having reversed polarity, whilst 
the anti-hypercomplex component corresponding to X2 reversing sign is projected/converted into a hypercomplex 
signal: thus the three anti-hypercomplex components are not individually annihilated like in the 1-D case. Two of 
the anti-components thus cancel, whilst the third is projected into a hypercomplex object. These results will of 
course be necessary to derive the properties of the hypercomplex wavelet transform of the phase-shifted signal. 



E. 9-Hypercomplex Extension of Phase-Shifted Signal 

Note that the 6'-hypercomplex extension of the phase shifted signal has the FT: 

.F{(Ai»^++^} = (1 +sgn([r_,f]2)) (1 - fcsgn([r_,f],)).F{(Ai»} (r_,f) 



(1 +sgn([r_,f]2)) (1 - ksgn{[r.gi],)):F{{Alf^,g)} (f). 
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The FT of the separable phase-shifted signal is 



/oo poo 
/ 
-oo J — 



oo J — oo 
-27rj(r_ex)^r_ef fi 



COS (27ra_5i(r_0x) — ^^^i) cos ( 27r/3_5i (r_6ix) — 



= ^ [e-^'^-Hl+sgn([r_ef]i))+e^'^-Hl-sgn([r_ef]i))] G_e(r--ef) 

e-^'''-^(l + sgn([r_ef]2)) + e^'^^'^l " sgn ([r^^fls)) ■ 
After some algebra we then obtain: 



{++)- 



(f) = (l + sgn([r_ef]2))(l-fcsgn([r_ef]i)) 



_ ^e-^-^=.i(l + sgn([r_ef]i))+e^'^-Hl-sgn([r_ef],))J G„e(r- 
e-^'''-^(l + sgn([r_ef]2)) + e^^-\l - sgn ([r^efJz)) 

Hence we may deduce that, as would be expected for g{-x) a separable signal we determine that: 



(++) 



(x)=e-^^-^i++;(x)e-^-^= 



(A-28) 



F. Hypercomplex CWT of a Phase-Shifted Signal 

We note that the hypercomplex CWT of the observed signal g((x) is by equations (IA-26t and (IA-28t : 



w 



{++) 



w 



(ee) , 



(++) 



-30s 



from equation (IA-28t . Thus the result follows. 



G. FT and QFT of the Monogenic Wavelet Transform 

Recalling equations ( I37t . (1^^ as well as (1^^ we obtain the FT 

wii^\C.g) = <^(C;5)-i<^(C;5)-jwf (C5) 



[1 — k cos I 



+ sm{4>t-9)]W';>{C9)- 



(A-29) 



The corresponding QFT is 

/(+) 



+ 



i-fc 

2 

1 + fc 
2 



[1 + cos (i/fc - 0) + sin (i/fc - 0)] G(q6)a^('-)* (ar^ecib) 



1 + cos (i/fc + e)+ sin (i/fc + 0)] G( J^/2qfe)a^'('^)* (ar-„e^^/2q6) 
^)'q('^; 5) + cos e [cos i/fe + sin Ub] 5) + sin 9 [cos t't - sin u^] iW^''^{K; g)j. 



(A-30) 
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H. Wavelet transform of 9-Monogenic & Anti-Monogenic Decomposition Components 



We find that the Fourier transform of g^^l 0,(x) is given by: 



{+) 



(1) 

-e',e',R- 



-(2) 



[l + sin(0{, — 9) + cos(0{, — 9) cos(0{, — 9') + sin(0f, — 9') + sin((/);, — 6) sin{(j)i, — 9') 



+ A; (— cos((/){, — 9) — cos( 



sm 



I cos 



We remark that with 6 = 9' this then becomes 



(+) 



{2 + 2sin(0b 
{2 + 2sin(0b 



9) -2k 


cos( 




9) -2k 


cos( 




9) -2k 


cos( 





+ cos(0b - 0)sin(</), - ^0)] W^^\c,g). 

(A-31) 



^)} <^(C;5) 

^)}G(f,)a^W*(ar_ef,) 



= 2G%(f,)a^W*(ar_,f,) = 2W^('-)(C;5i+,;, 
Furthermore mutatis mutandis with g^Sgg{x.) replacing g^^Jg{x.) we derive that 



(+)^ 



0. 



(A-32) 



(A-33) 



This completes the proof. We can note that equation ( IA-3U establishes the necessity of introducing the ^-monogenic 
signal. The Monogenic wavelet transform will not annihilate the monogenic signal for example, and in the subsequent 
section this may produce problems. 



/. Construction of Local 9-Monogenic Signal 

We note the form of the Fourier transform of w^\^]g) from equation ( 11061 ) and then from equation (H^ we 
find that (recalling 



[a,0,hf)- 



w, 



(r) 



(+) 



w 



(r) 



(+) 



[1 + sin 
[1 + sin 

i 



- 9) - fccos(4 - 9)] G^e (r-efb) a^^''> (ar^efb) 
i^b - 0) -kcosi^b - 0)] G(ffe)a^M* (ar_effe) 



Thus the monogenic wavelet transform coefficients correspond to the monogenic extension of the real wavelet 
transform in the rotated frame of reference. For the second equality we start by using equation ( H3t and take 

5(x) = I 



X, 



By the linearity of the wavelet transform we have that with 9' = 9 using equations (IA-32t and (IA-33t : 



1 

2 L 
1 



(A-34) 



and thus the result follows. Hence the operation of constructing the monogenic wavelet transform can either be 
considered in light of the operation of 1) taking the monogenic wavelet transform of a real signal, 2) finding the 
real wavelet transform of a ^-monogenic signal or as finding the 6'-monogenic image of the wavelet transform in 
b of the rotated by 9 signal. 
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J. Local Invariance under Rotation 

If the monogenic wavelet is constructed from an isotropic wavelet, then it follows that: 



(1) 



/oo poo roo /"OO 

/ <)(C;9)e^^'-^^^^d2f.= / / G(f,)aM/ (ar_,f,) e^^-^f^ d^f, 
-OO J —CO J —OQ J —OO 

/OO /"OO 
-OO J —OO 



(A-35) 



where we have used that the real mother wavelet is isotropic in (1). In a similar fashion, we calculate the for the 
CWTs with respect to the first and second Riesz components of the monogenic wavelet. 



cos(^)w;J,') (^o; 5) + MO)wf (^o; 5), 
- sin(0)u;J) (6; 5) + cos(0)u;Jf) (6; 5)- 



Using the above expressions (IA-35t . and (IA-36t . trivially follows that 



wl' ' {i\ g) 



w, 



(2)2, 



W, 



(+) 



(6; 5) 



(A-36) 



(A-37) 



Hence, if an isotropic real mother wavelet is used, the magnitude of the monogenic wavelet coefficients is invariant 
to to the choice of orientation. 



K. Directional Selectivity 

For simplicity we model the signal is directional in the Fourier domain, i.e. via equation jl 141) . Note that the 
wavelet transform corresponds to patterns at period /. It then transpires that: 

f OO fOO 

Tir(l)/> X 2i7rfJb ,2, 



/OO roo 
-OO J ~OD 



— OO J —CO 

Jo Jo ^ 

roo ^ 

= cos{9-(j)o) / fbG{fb)a¥^\af) sin {2tt fbb cosixb - M) dfb, 
Jo 

In a similar fashion, 

/•OO ^ 

it g) = sin(0 - 0o) / /6G(A)aM/(^)(aA) sin (2^A6cos(xb - 4)) 4fb- (A-38) 
Jo 

For (1) to hold, V^'^H^) = V'^'^H^)' isotropic. For a directional signal, when w^^^ (^;5) 7^ 0, 

^ - 4 = tan-i f ^1^1 . (A-39) 

Thus the local directionality of a directional signal, can at any rotational angle of the wavelet transform, be 
determined from the isotropic monogenic wavelet transform. Furthermore it is clear that 

w^!U$^;g) =0, (A-40) 



and hence by equation (11071) . the stated result follows. 
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L. CWT of Phase-Shifted Signal 

For br 
note that 



For brevity we denote the wavelet transform of the phase-shifted signals by w^^j i^]9) = Ag^g). We 



Ae.5(x) 



cos (27r(/)(x) — 9s 



1 r 
2 



-Bi „Jx)0, ~(+) / \ , ,Jx)e, ~(~) 



Hence with the additional assumption that ej;_|,g(x) = ep^g(b) is constant over the spatial width the wavelet 
transform is averaging over, i.e. |x — b| < axi we retrieve that 

1 



e---«,«(^)^^4+)(^; (X)) +e--.^(^)^^4+)(^;5^W (X)) 

where we used equation (IA-33t . Also from equation (IA-32t we may note 



2 
1 



(A-41) 



(A-42) 



and thus the result follows. 



